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1. Extraordinary Astronomy 


There are few things more intriguing in high energy astrophysics than the study of the highest 
energy particles in the universe. Where and how these particles achieve their extreme energies 
is of interest not only to the astrophysicist but also to the particle physicist. At GeV and TeV 
energies the problem is manageable since the physics is known and the acceleration processes 
feasible. But the energy spectrum extends to 10 20 Ev and there the problem of their origin is 
both more difficult and interesting; in fact at these high energies we do not even know what 
the particles are. 

The study of the origin and distribution of relativistic particles in the universe has been a 
challenge for more than 80 years but it is only in recent years that the technology has become 
available to really address the question. Although something can be learnt from studies of 
composition and energy spectrum, the origins (and thence the acceleration mechanisms) can 
only come from the direct study of the neutral particle component (in this respect the high- 
est energy particles are effectively neutral since they are virtually undeflected). The feasible 
channels of investigation are therefore the study of the arrival directions of: (1) TeV photons 
(covered by the following U.S. experiments: STACEE, Whipple/ VERITAS, MILAGRO and, to 
some extent, by EGRET/GLAST); (2) neutrinos of TeV energy and above (AMANDA/KM3); 
(3) the highest energy cosmic rays (HiRes, Auger). 

While these studies represent a form of astronomy they are the astronomy of the extraordinary 
universe , the universe populated by the most dynamic and physically exciting objects, the 
universe of the high energy astrophysicist whose cosmic laboratories represent conditions beyond 
anything that can be duplicated in a terrestrial laboratory. This extraordinary astronomy may 
say little about the normal evolution of stars and galaxies but it opens windows into cosmic 
particle acceleration where new and strange physical processes take place. 


2. Role of Whipple Collaboration 


The Whipple Gamma Ray Collaboration was formed in 1982 with the aim of improving the 
sensitivity of the ground-based atmospheric Cherenkov technique by the use of imaging cameras 
on large optical detectors. The approach was successful and the discipline of ground-based TeV 
astronomy was established. The first confirmed galactic source was discovered by the Whipple 
Collaboration as were the first two extragalactic sources. The technique has now matured and 
several new and important astrophysical results have been produced. 

The 7-ray group of the Smithsonian Astrophysical Observatory is a key component of the 
Whipple Gamma Ray Collaboration. With all its members resident in Arizona, all its activities 
are devoted to research in ground-based 7-ray astronomy. 

3. The Crab Nebula 


TeV observations of the Crab Nebula made at the Whipple Observatory over the last decade 
indicate no evidence for flux variability in the energy range from 400 GeV to more than 10 
TeV, making this object a valuable standard candle for calibration purposes. 



The spectrum of the Crab Nebula exhibits a remarkably broad dynamic range, spanning the 
energy range from less than 10 -4 eV to at least 5xl0 13 eV. Gould (Gould 1965) postulated that 
the entire spectrum (then not even known to extend to X-ray energies) could be explained 
by a Compton-synchrotron model. From the lowest radio frequencies up to sub-GeV energies, 
the spectrum is dominated by the synchrotron emission of relativistic electrons with energies 
extending up to 1 PeV; these are now believed to be accelerated by the pulsar wind. The 
GeV-TeV photons arise from inverse Compton scattering of synchrotron photons or other low 
energy photons (e.g. the microwave background) in the nebula by the relativistic electrons. 
Over the dynamic range 500 GeV to 5 TeV the spectrum is well represented by a power law 
of the form J = (3.2 ± 0.7).10 -7 x (£/lTeV) 2 ' 49±0 ' 06±005 m -2 s -1 TeV -1 (Mohanty et al. 1998). 
A recent observation by CANGAROO at low elevations provides data in the energy region > 
7 TeV, which suggests that a power law spectrum might extend up to 50 TeV. However, the 
air shower array limits at energies greater than 100 TeV fall below' the extrapolation. If this 
linear power law is extrapolated to lower energies, it passes more than a decade above the upper 
EGRET point at 5 GeV. A form which is quadratic in Log (E) (Hillas et al. 1998) satisfies 
both the GeV and TeV data and is consistent with upper limits at energies greater than 100 
TeV. This spectrum implies a magnetic field in the nebula of 160/zG, close to the equipartition 
value. 

4. Supernova Remnants: Shell-type 


The canonical theory of cosmic ray origins suggests that they emanate in shell-type SNRs. 
High energy 7-ray observations can indicate which SNRs have a large content of relativistic 
cosmic ray hadrons, and TeV 7-ray observations, in particular, have the sensitivity and angular 
resolution to reduce background confusion. EGRET measurements of SNRs are not definitive 
because the detector has low angular resolution at 100 MeV and measurements are masked by 
gas clouds. Collisions of cosmic ray nuclei with the interstellar medium result in the production 
of neutral pions which subsequently decay into 7-rays. These processes result in a secondary 
7-ray spectrum which follows the primary cosmic ray spectrum at energies above ~10 GeV up 
to ~1/10 of the maximum proton energy of ~ 100 TeV. 

Calculations indicate that the luminosity of nearby SNRs should be sufficient for detection by 
the most sensitive VHE 7-ray telescopes. If there is a density enhancement from a molecular 
cloud, current ACTs and EGRET should already be able to detect the 7-ray emission from 
some objects. For a source spectral index of a=2.1, Drury et al. (1994) estimate the integral 
7-ray flux at earth to be 


F(> £) « 9 x 10 


-11 


(1 


E 


TeV 


-1.1 


/ _0Esn_\ f 
\10 51 erg/ Vl 


d 


kpc 


)"( n ==0 


cm 2 s 1 


(1) 


where 9 is defined as the fraction of the supernova energy Esn converted into cosmic rays, d is 
the distance to the SNR, and n is the average density of the ISM around the remnant. 
Assuming 9 « 0.15, a fairly conservative value for the average density of n « 0.2cm -3 and 
the canonical value E S n ~ 10 51 ergs, gives a flux F{> 200 GeV) « 1.6x 10 -11 (d/l kpc) -2 - 
photons cm -2 s -1 . This flux lies close to the sensitivity limit of the current generation of TeV 
telescopes. 

For several of the unidentified EGRET sources there is now' evidence for an association with 
a SNR (7-Cygni, IC443 and W44. Observations by the Whipple telescope do not confirm the 



Table 1 : SNR Upper Limits 


Source 

EGRET Flux 
xl0 _7 cm~ 2 s _1 

Whipple Upper Limit 
xl0~ n cm~ 2 s _1 

W44 

5.0 

<3.0 

W51 

<3.2 

<3.6 

7-Cygni 

12.6 

<2.2 

W63 

<1.9 

<6.4 

Tycho 


<0.8 

IC443 

5.0 

<2.1 


expected extension of the 7-ray spectrum[?]Buckley et al. 1998). Some of these limits are listed 
in Table 1. They lie a factor of ten below the predicted spectrum, but there is still significant 
uncertainty in the remnant parameters so that the canonical model is not yet in jeopardy. It 
is also possible that the EGRET associations with the SNR are incorrect. Although not yet 
definitive, the upper limits require that the source spectrum for IC443 and 7-Cygni be steeper 
than E ~ 2A ; this would imply that the additional steepening of the cosmic ray spectrum due to 
propagation effects in the galactic disk should be —0.3 rather than the preferred value of —0.6. 


5. Pulsars 


A likely source of TeV 7-rays is radio pulsars; this is particularly so since more than six of 
them have been detected as 7-ray emitters by CGRO. There were also unconfirmed reports of 
pulsed emission from several early TeV experiments. The most convincing detection, a steady 
pulsed signal from the Crab pulsar, was from the Durham group (Dowthwaite et al. 1984); the 
observed flux was ~ 1 x 10 -11 cm -2 s -1 at 1 TeV. More sensitive observations show that less 
than 5% of the total TeV flux is pulsed (Gillanders et al. 1997). Evidence for the detection of 
a pulsed signal was also suggested for the Vela pulsar and for the Geminga pulsar. However, 
observations with more sensitive imaging ACTs failed to confirm these early results. 

An extensive search by the Whipple group for emission for the two other EGRET pulsars 
(Geminga: Gillanders et al. 1997; PSR 1951+32 Srinivasan et al. 1997a) that are visible in the 
northern hemisphere failed to find any evidence for emission. No TeV pulsed signal has been 
detected from the other (southern) CGRO pulsars, PSR B1706-44 and PSR B 1055—52 either. 
No TeV observations have been reported yet on PSR B0656+14. 

The upper limits for the TeV flux from the EGRET 7-ray pulsars are below the fluxes extrap- 
olated from the GeV region of EGRET detection assuming a power law spectrum of constant 
index. The emission from the pulsar magnetosphere seems to turn off or to fall off steeply in 
VHE region. It has been have argued that the pulsed component will extend to VHE energy 
in outer gap models. 

The analysis of pulsars is not trivial since it is necessary to ascertain that the system timekeeping 
is sound, that the absolute timing is accurate, that the pulsar emphemeris (usually based on 
radio observations) is correct, and that the times of arrival are barycentered and folded correctly. 



Table 2: Whipple Limits from Gamma Ray Pulsars 


pulsar 


EGRET HE flux 
(10~ 7 cm -2 s _1 ) 

Group 

VHE Flux 
(10 -12 cm -2 s -1 ) 

Eth 

(TeV) 

Crab 

unpulsed 

7.7 ±0.8 

Whipple 

8.8 (E/TeV)" 1 - 69 

0.4 


pulsed 

23 

Whipple 

< 0.2 

0.25 

Geminga 

unpulsed 

- 

Whipple 

< 8.9 

0.5 




Whipple 

< 5 

0.5 

1951+32 

pulsed 

1.6 ±0.2 

Whipple 

<5.4 

0.3 


Once a periodic 7-ray pulsar signal has been detected, then all systems are verified; however a 
null result leaves the suspicion that one or more of the above may have been incorrect. 

To verify all components of the Whipple pulsar data taking and analysis, the telescope was 
converted into an optical telescope to observe the optical pulsar in the Crab Nebula. This 
was achieved by stopping down the central pixel of the 10m camera and recording the output 
on a scaler readout every 1 millisec. Folding the data taken over a ten minute observation at 
the predicted pulsar period revealed the familiar double-peaked light- curve of the Crab pulsar 
with the main peak exactly at the predicted phase (Srinivasan et al. 1997b). The null results 
reported above can therefore be treated with confidence. 

6. Unidentified EGRET Sources 


Despite extensive searches for counterparts at long wavelengths for the 30 odd low latitude 
unidentified EGRET sources close to the galactic plane, no identifications have been made and 
their nature is unknown. A search for counterparts at high energies is justified because many of 
the sources exhibit flat spectra and the imaging ACTs have good source location capabilities. 
The Whipple telescope has been used to observe a number of these (Buckley et al. 1997) but 
no significant emission has been seen. 

7. Bursts 


A feature of 7-ray astronomy has been that as the energy is increased there is an increase in 
the degree of temporal variation. As seen above, in AGN studies, TeV variations with doubling 
times as short as 15 minutes have been observed. It would not be unexpected that TeV 7-ray 
bursts would be observed, either as the tail end of classical 7-ray bursts or as a manifestation 
of a new phenomenon. 

There are a number of exotic suggestions that justify a search for 7-ray bursts at TeV energies. 
These include emission from the decay of primordial black holes(PBHs) and cosmic strings. 
Most models predict a final explosion of energy when all possible evaporation channels are 
available, but the number of degrees of freedom of emission is a highly controversial issue. 
Upper limits have already been set to the PBH density by atmospheric Cerenkov and air 
shower array experiments under various assumptions. 



Predictions have been made that (i) non-conducting cosmic strings acquire cusps that are 
smoothed out by emitting bursts of TeV 7 -rays over ill-defined time-scales and that (ii) su- 
perconducting strings with a saturated current produce a jet of fermions which decay to TeV 
7 -rays over a 1 second period. In general these predictions have sufficient free parameters that 
non-detections are not serious limitations; however the exciting new physics that a positive 
detection would indicate fully justify the search for new burst phenomena as TeV detectors are 
developed with improved flux sensitivity. 

Since August 1992, the Whipple collaboration has searched for bursts on a one second time- 
scale, both on-line and in the archival data-base. The methodology required for a search of this 
nature (where the source can come from anywhere in the field and the onset is unknown) is 
unlike that used in the search for point sources and required the development of new analysis 
techniques (Connaughton et al. 1998). The number of bursts of 3 or more per 1 second was 
compared with the expectation value and no significant excess was obtained over the 4 years of 
the data-base. From this null result a reliable limit to the density of PBHs could be derived. 
Although the serendipitous overlap of an atmospheric Cerenkov telescope with its limited field 
of view with a classical BATSE- detected burst is unlikely, there is some hope of detecting the 
delayed high energy component seen in some bursts. It is possible to make rapid follow-up 
observations of BATSE bursts using source positions distributed on the BACODINE network. 
The effectiveness of these observations is limited by the restricted duty-cycle of ACTs, by 
the imprecise position locations of the bursts, by the slew speed of the telescope and by its 
limited field of view. Nonetheless observations by the Whipple Collaboration of 16 BATSE 
positions, one acquired within two minutes of the reported BATSE burst time, have been 
reported (Connaughton et al. 1997; Boyle et al. 1997). However in no case did the FOV of 
the telescope overlap the complete error box in source position uncertainty. No evidence of 
TeV emission is found and upper limits to the high-energy delayed or extended emission were 
derived based on assumptions of the source positions. 
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Abstract 

Although atmospheric Cherenkov telescopes have restricted fields of 
view their fluence sensitivity warrants a search for gamma-ray burst phe- 
nomena A search for 400 GeV gamma-ray bursts on a Is time-scale 
using archival data taken between 1988 and 1992 with the Whipple Col- 
laboration 10m reflector is presented. No evidence of such bursts is found. 
Bursts of TeV gamma rays have been predicted from exotic astrophysical 
objects such as Primordial Black Holes and Cosmic Strings. An upper 
limit to the number density of exploding PBH of 3.0 i 1.0 x 10 6 pc^yr” 1 
is calculated. 

Contact: Valerie Conn aught on , ES 84, Marshall Space Flight Center, AL 
35812. vc@msfc.nasa.gov. tel: (205) 544 6494 FAX: (205) 544 5800. 
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1 Introduction 


The detection of gamma-ray bursts with peak luminosity at MeV energies has 
been one of the most exciting discoveries in high energy astrophysics [9]. Al- 
though still unexplained after 30 years of intense research the phenomenon has 
inspired searches for counterparts at a variety of wavelengths [14]. It has also 
opened the possibility that there might be similar phenomena that have peak 
luminosity at other wavelengths; here we consider a TeV search for such phe- 
nomena on time-scales of one second using the Imaging Atmospheric Cherenkov 
Technique. 

The Whipple Observatory Gamma Ray Collaboration has been involved in 
searches for gamma-ray bursts since 1978; initially multi-element first genera- 
tion systems over long baselines were used [17]. The development of the Imaging 
Atmospheric Cherenkov Technique has significantly improved the sensitivity for 
the detection of point sources [25]. The Whipple 10m imaging camera has 
successfully detected TeV emission from several point sources - the Crab Neb- 
ula [24], Markarian 421 [19] and Maxkarian 501 [20]. It has been the most 
sensitive instrument of its kind for source-centered observations but its off-axis 
sensitivity has not so far been fully characterized [1]. 

Here the response of the Whipple 10m reflector to non-source-centered gamma 
rays is assessed using shower simulations, and a method is developed to search 
for gamma-ray bursts of one second duration. We report on a search through 4 
years of the Whipple database. 

There are a number of exotic suggestions that justify a search for gamma-ray 
bursts in this parameter space. These include emission from the decay of pri- 
mordial black holes (PBHs) [13] and cosmic strings [22], The cosmological and 
physical importance of PBHs is well-established [12]. As the PBH evaporates, 
its temperature increases in proportion to the mass loss. Most models predict 
a final explosion of energy when all possible evaporation channels are available, 
but the number of degrees of freedom of emission is a highly controversial issue, 
and ranges from the conservative standard elementary particle model, where it 
reaches a maximum once the three generations of fermions are free to be pro- 
duced, to the runaway Hagedom model where the number of emission modes 
increases exponentially beyond this [12]. For the conservative model an opti- 
mum time-scale of 1 s is found for observing the explosion of a PBH using the 
Whipple 10m telescope. During this final second 9.1 x 10 28 photons are emitted 
above 0.4 TeV. 

The fluence sensitivity of the detector over this time-scale is of order 10 -9 
erg cm -2 which is comparable to or better than that achieved in most bands 
of the gamma-ray spectrum. The sensitivity of the technique is compared to 
that achieved in searches by other gamma-ray experiments, operating at energy 
thresholds higher than the Whipple instrument. 
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2 The Whipple 10m Telescope 

Situated at an altitude of 2.3 km on Mt.Hopkins in Arizona, the telescope 
operated by the Whipple Collaboration has been used as an imaging device 
since 1982; the camera has been improved from a 37 to a 109 element imag- 
ing system during that time. The telescope consists of a 10-metre dish (with 
248 front-coated mirrors) on an alt-azimuth mount and a camera containing 
109 photomultipliers (PMTs). The inner 91 PMTs are 1.1cm diameter tubes 
(Hamamatsu R1398) and are connected via amplifiers to trigger discriminators. 
An outer ring of 18 tubes is independent of the trigger logic and contributes 
only to the imaging of the event. Until 1993, these outer tubes were 5cm tubes 
giv ing a total field-of-view of 3.75°. They have been replaced by the smaller 
tubes so that the camera field-of-view is now 3°. The arrangement of PMTs in 
the focus box prior to 1993 is shown in Figure 1 and a detailed description of 
the instrument can be found in [5]. 

The cam era is triggered if the signal in two or more PMTs exceeds a preset 
threshold in a 10 nsec time interval. When the system is triggered i.e. an event 
is registered, the light level in each PMT is recorded over a 25 nsec gate width. 
The trigger rate during the epoch of the observations reported here was typically 
8 events per second. The time is derived from a GPS dock and is recorded with 
event with an accuracy of 0.1 msec. The deadtime is 1.5 msec and the 
distribution of time differences of successive events from a 60 minute data file 
taken at a telescope elevation of 80° is shown in Figure 2. 

The digitized images of the Cherenkov light are analyzed off-line. A typical 
image recorded by the present camera is shown at various stages of processing 
in Figure 3. Moment-fitting routines are used to define each image by an ellipse 
whose parameters define the event (shown in Figure 4). The event is accepted as 
a potential gamma-ray event or rejected as a hadron-initiated event depending 
on the values of these parameters. The axis of the image is defined by minimizing 
the signal-weighted sum of the squares of the perpendicular angular distances of 
the pixels. The rms spread of light perpendicular to and along this axis define 
the semi- min or and semi-major axes of the ellipse and these are known as the 
width and length of the image. Image compactness is defined by concentration, 
which is the fraction of the sum over all non-zero pixels of the cleaned image 
in digital counts (known as size ) contained in the two or three brightest pixels. 
The other parameters shown in Figure 4 relate to the orientation of the image 
axis relative to the source position in the field-of-view. 

If the position of the source is known (usually at the centre of the field of 
view) then the orientation of the image relative to that direction is important 
for source detection. A combination of shape and orientation parameters has 
been used by the Whipple Collaboration to reject 99.7% of recorded background 
while keeping 50% of gamma-rays; the set of selection parameters that are most 
useful for observations of a point-source in the centre of the field of view are 
called Supercuts and are defined elsewhere [21]. A comparison of the width and 
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length d omains for photons (simulated) and real background (data taken with 
the telescope pointed at the zenith) can be seen in Figures 5. 

The success of Supercuts and the Imaging Atmospheric Cherenkov Technique 
in the detection of point sources owes much to Monte Carlo simulations of 
photon and background hadron-initiated cascades. Owing to the impossibility of 
testing the atmospheric Cherenkov detector at an accelerator, analytical models 
and simulation of phenomena are used to optimize instrumental design and in 
the development of analysis techniques. Data taken during observations of the 
Crab Nebula, which plays the role of a TeV gamma-ray standard candle, are 
used to perfect new techniques. 


3 Off-axis Sensitivity of 10m reflector 

The characteristics of the 10m reflector in point-source mode of operation are 
well understood and defined [5]. What is less well characterized is the sensitivity 
of the instrument when an area of the sky is probed for a source without know- 
ing its exact position, i.e., where each point in the field-of-view is a potential 
source. It has been demonstrated experimentally that the system is sensitive to 
a point source of gamma rays up to 1° off-axis [1], [10]. In the work reported 
here the response over the full field-of-view of the 10m reflector is investigated, 
linin g simulations, in order to use the instrument as a relatively wide field burst 
detector. 

3.1 The off-axis Triggering and Imaging of Gamma Rays 

The geometrical field-of-view of the 10m camera between 1988 and 1992 was 
3.75°, but the imaging technique allows reconstruction of images whose origins 
lie outside this area. Conversely, some of the images of showers from sources on 
the edge of the camera face may fall outside the field-of-view so that the camera 
has reduced collection area for such sources. 

In order to assess the sensitivity of the instrument for sources lying at var- 
ious distances from the centre of the camera, a database of nearly 8000 sim- 
ulated gamma-ray- induced showers of energies ranging from 0.2 TeV to 0.8 
TeV was compiled using a simulation program based on the KASCADE code 
of Kertzmann and Sembroski [15]. The Whipple collaboration has used several 
different simulation codes in the various research institutions (e.g. MOCCA at 
the University of Leeds) and comparison of the simulations shows that they are 
in agreement with each other to within 10% in the Cherenkov light yields and 
image parameter distributions that they produce. 

The response of the instrument to the simulated gamma-ray showers is 
gauged by mapping each Cherenkov photon in a shower from the mirrors onto 
the imaging camera and building up an event for each shower. 
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Displacing the focus box from its position at the focus of the reflector and 
constructing events as they appear to this displaced camera is equivalent to 
simulating the response of the camera to sources at an offset from the centre 
equal to the displacement of the camera but in the opposite direction. The 
response of the camera to gamma-ray showers from sources with 0° to 2.5° 
offsets from the centre was investigated by successively displacing the camera 
at 0.25° intervals along an arbitrary line, 17° from an axis through the centre of 
a row of PMTs. A different source displacement angle may be chosen without 
affecting the results. The treatment of simulated images is identical to that 
subsequently applied to real data. The parameterization and cutting procedures 
execute rapidly relative to the shower and instrument simulations, and it was 
possible to experiment with a variety of parameters and analysis methods before 
deciding on a suitable gamma-ray burst search method. 

As the source position moves further away from the centre of the camera, 
less Cherenkov light falls on the detector at a given impact parameter. The 
software trigger is a requirement that two of the inner 91 tubes register at least 
40 photo-electrons and the percentage of events triggering the system at each 
impact parameter is shown in Figure 6 for source offsets of (a) 0°, (b) 0.75° and 
(c) 1.5°. 

By weighting the fraction according to impact parameter value (the collec- 
tion area is bigger at larger impact radii), a value for the trigger efficiency over 
the field of view, equivalent to the effective area of the detector for each source 
offset, is calculated: 

Ei ZvrifaSr 

Trigger efficiency = - 1Q4 — (1) 

where r< (=i) is the radius of the annulus at impact parameter i, thickness 
Sr = 2 metres, and <j>i is the fraction of the 20 showers at radius r, which 
would trigger the telescope. The summing integer i is incremented from 12 to 
200 in steps of 2. At impact parameters below 10 metres, the showers tend to 
saturate the central tube, and the Cherenkov light pool falls off rapidly beyond 
150 metres, so that we define an ideal detector as one which is triggered by all 
showers with impact parameters between 10 and 150 metres, and which fails to 
register any events falling outside that area. The trigger efficiency is divided by 
7 x 10 4 m 2 to give the efficiency relative to this arbitrary ideal detector. An 
efficiency of greater than 1 is, therefore, possible since showers lying beyond 
150 metres may trigger the system. This procedure was implemented for all 4 
shower energies, and the results are shown in Figure 7. Above threshold (> 0.4 
TeV) the instrument would register over half of the gamma-ray showers from 
sources out to 1° away from the centre, but it is obviously more efficient at 
detecting those nearer the centre than at large offsets. 

Image parameters have traditionally been used, in some combination of 
shape and orientation cuts, to discriminate against the hadronic background. 
In this search, however, only the shape parameters are relevant since there is 
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no preferred direction due to the unknown position of the source in the field-of- 
view. The most successful shape selection was found to involve cuts in length 
and width , using the values derived for Supercuts: 

0.073° < width < 0.15° 

0.16° < length < 0.30° 

Figure 8 shows the percentage of 0.4 TeV events at offsets (a) 0°, (b) 0.75° 
and (c) 1.5° which trigger the camera and satisfy these image shape require- 
ments. Combining the percentages for each impact parameter produces the 
results which are shown in Figure 9. Above the detector energy threshold, the 
collection area for gamma rays is fairly flat, and showers from sources which 
are up to half a degree off-axis are still selected over 50% of the time. There 
are two competing effects which explain the shape of the curves: at high en- 
ergies, showers with large impact parameters produce enough Cherenkov light 
to trigger the telescope more often than distant low energy showers, but the 
Supercuts selection is biased towards small showers, so that the higher energy 
showers are less likely than the low energy events to survive the width and 
length cuts regardless of impact parameter. At large source offsets the fraction 
of events selected decreases with shower energy even though the high energy 
events appear to trigger the camera more often than at lower energies. If the 
showers triggering the system at these large offsets are examined, it can be seen 
that they are the longest images in the high-energy databases. 

Observations of the Crab Nebula at angular offsets from the centre of the 
field-of-view of the 10m telescope were used to confirm the efficiency curves 
derived from these simulations. The Crab Nebula is a steady source of TeV 
gamma rays, and by comparing the rate of gamma rays measured after applying 
the Supercuts selection technique to data from September 1996 at offsets ranging 
from 0° to 1.5°, one obtains the efficiency for gamma ray collection displayed 
in Figure 10. The errors in Figure 10 reflect the variation in event rates for the 
offset observations, and are large owing to the small exposure time: only 1 hour 
of source observations were made at a source offset of 0.5° from the centre of 
the field-of-view, 1.5 hours at 1.0°, and 2 hours at 1.5° offset. Because the field- 
of-view of the telescope was smaller when these offset observations were made 
compared to that simulated in this study, a perfect correlation is not expected, 
and one anticipates that the efficiency at larger source offsets be slightly lower 
relative to source-centered observations than in the simulations. 

3.2 Orientation of gamma-ray images 

The orientation of the ellipse fitted to each image is represented by its major 
axis, and the most likely point-of-origin of the shower progenitor on the field- 
of-view lies on this axis at a distance d in degrees related to the ellipticity of 
the image: 
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d = 2 — 2(width/ length) (2) 

This algorithm was developed based on simulations of gamma-ray showers 
from sources at the centre of the field-of-view, and was found to be accurate 
to about 0.3° either side of this point [1] for source-centered observations. The 
relation between the image, its centre, and the most likely point-of-origin are 
shown in Figure 11. 

Because the ellipses derived from the moment-fitting routines are symmet- 
rical, the point-of-origin may lie either side of the centre of the image. The 
following equations are used: 


vW - Veen) 2 + (x 0 r - x cen ) 2 = 2 - 2 x width/length 

Vor = mXXor+C 

to find the most likely points-of-origin i OT , y or on both sides of the centre of 
light x cen , y C en along the major axis y = m x x + c. The distance 


\/(i/0 - Vor)' 1 + (zo - lor) 2 (3) 

where xq, yo are the true positions of the simulated source, is a measure of how 
well Equation 2 finds the point-of-origin of a shower. The tolerance is defined 
as the maximum allowed value of this distance, in degrees. A range of tolerance 
levels from 0.1° to 0.7° was explored to determine the off-axis efficiency of the 
point-of-origin determination, and the results are displayed in Figure 12. The 
total efficiency $ for each energy was obtained by combining source positions 
and weighting each efficiency (4> r ) according to source offset (r). The efficiency 
values featured are relative to an ideal detector with 4> r = 1 at each source 
offset where the gamma-ray images are selected and their origin is correctly 
determined so that $ = '52<t> r r / with r varying from 0.25° to 2.5°. The 
efficiency curve is fairly constant as the source is moved away from the centre 
and appears independent of impact parameter and photon energy. 

Real background data were subjected to the same analysis so that the effect 
of varying the point-of-origin tolerance value on background acceptance could 
be assessed. It was found that in looking for bursts of 3 or more events over 
a 1-second time-scale, the highest signal to noise ratio could be obtained with 
a tolerance level of 0.425°. This selection criterion enables 66 ± 5% of shape- 
selected simulated events (over the impact parameter range 12 to 200 metres and 
source offset 0° to 2.5°) to be correctly located. The shaded areas in Figure 11 
show the regions from which the imaged events could have originated with this 
tolerance value. 
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It is then possible to apply this calculated off-axis response of the detector to 
the search for counterparts to the delayed component of B ATSE bursts (reported 
elsewhere [7]) and to the serendipitous detection of gamma-ray bursts during 
the normal operation of the'gamma-ray telescope (described below). 


4 Searches for 1 s TeV Gamma-Ray Bursts 

The archival data of the Whipple Observatory taken with the 10m reflector 
with an energy threshold above 0.4 TeV between September, 1988 and Septem- 
ber, 1992 were used for this search. In general these observations were in the 
form of files of (approximately) 28 minute duration containing information on 
each event (the arrival time, the outputs of the 109 ADCs from the 109 pixels, 
and other housekeeping data). Examinination of the observing logs allowed the 
identification and rejection of 1% of the data including (i) observations made 
when there were instrumentation problems (ii) data taken with discriminator 
thresholds or high-voltage values significantly different from those used in nor- 
mal operation, and (iii) observations made while testing new or experimental 
configurations (for example, experiments with filters). A further 1% of the files 
were rejected because of obvious anomalies, such as a large number of bright 
stars in the field-of-view, or because they were taken at a telescope elevation less 
than 35° above the horizon. Following rejection of unsuitable data, 2217 hours 
of observations comprised the database. The database was divided according to 
the event rate of candidate gamma rays. Variation in event rates was caused by 
differences in zenith angle, mirror reflectivity and discriminator settings. The 
elevation of the source under observation is known to affect the energy threshold 
of the detector, and also influences the appearance of both gamma and hadronic 
images. A larger depth of atmosphere must be penetrated by the extensive air 
shower particles if they are to produce light lower in the atmosphere, so that a 
telescope pointed close to the zenith is sensitive to lower energy showers than 
if it is operated nearer the horizon. Consequently, the raw event rates decrease 
with increasing zenith angle - this relationship is shown in Figure 13 (top). De- 
spite this overall fall in event rates at low elevations, the Cherenkov image on 
the camera face is narrower on the horizon than at the zenith, and this effect 
results in more hadronic events being selected as gamma rays on the basis of 
the width parameter. Figure 13(bottom) shows the variation of image-selected 
event rate with telescope elevation. 

Nearly 80% of observations were made with the telescope elevation above 
55°, where data rates and image characteristics are fairly homogeneous. These 
observations form the bulk of the 1759 hours of data which had an average rate 
of < 0.2 candidate gamma-ray events per second, and comprise the largest of 
the three subsets of the Whipple database. The rest of the data was divided 
according to the average event rate; rate = 0.2 to 0.4 per second (394 hours) 
and rate > 0.4 events per second (64 hours). Standard routines were used to 
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flat-field the data, tubes with bright stars in their field-of-view were turned off in 
hardware or software, and the image parameters were calculated as described in 
[21]. The events were then characterized by these image parameters (and their 
arrival times). The evenness of the camera response across the field-of-view is 
illustrated in Figure 14, showing the distribution of event centroids in camera 
coordinates (degrees from centre) for a 60 minute observation. 

4.1 The Burst Search 

Initially the events were culled to include only those events whose images satis- 
fied the gamma-ray event criteria. Events which survived the width and length 
cuts described above were included in a reduced delta set which was submitted 
to a search for bursts of 3 or more events in a Is time interval. 

The arrival tim e of the first event was taken as the possible start of a Is burst 
and was compared to the time of the second event. If the time difference was 
under Is, the next event was examined. If not, the start of the time window was 
moved to the arrival time of the second event. The window was moved through 
the file until all events in a file have been processed, and the totcil number of 
bursts with from 3 to 10 shape-selected events in the file is obtained. No candi- 
date bursts were found with more than 4 events. The number of 3- and 4-fold 
burst events are shown in Table 1 in the column headed “obs”. A background 
(control) file was generated from each real data file by scrambling the events in 
each parameterized file, while maintaining the original time sequence. When the 
background file was subjected to the burst-search program, the image-selected 
events were the same as those picked from the real data file, but in a random 
order and attached to the original arrival times of other events. The number of 
3- and 4-fold bursts thus found is shown in Table 1 headed “exp”. The difference 
between the observed and the expectation is expressed in terms of the standard 
deviation, a defined as (obs + exp) 1/2 . In no case is the difference statistically 
significant. 


Table 1: Frequency of gamma-ray burst candidates in Whipple database vs. 


\JBBM 

Hours of 
Data 

3-event burst 

4-event burst 

exp 

obs 


exp 

obs 

a 

0.0-0.2 

1759 

4747 

4762 


234 

245 

+0.50 

0.2-0.4 

394 

13463 

13679 


1308 

1286 

-0.43 

> 0.4 

64 



-0.40 

1800 

1724 

-1.28 

Total 

2217 


28980 

+0.72 

3342 

3255 

-1.07 


To this point no use has been made of the fact that the events in the putative 
3- and 4-fold bursts must come from the same point in the sky. Hence all the 
bursts were subjected to the common origin selection described above. Less 
than 2% of the putative bursts passed this selection. These are listed as before 
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in Table 2; the expectation is derived by applying the common origin selection 
to the scrambled events. 

Again no excess of 3- or 4-fold bursts of candidate gamma-ray events on a 
1-second time-scale is found in the 2217 hours of the Whipple database. Two 
candidate four-fold events are shown in Figure 15. If a TeV component to a 
background of bursters exhibiting activity on a 1 s time-scale exists, then the 
number of these bursts is low and below the sensitivity of the instrument in its 
current form. 


Table 2: Frequency of gamma-ray burst candidates in Whipple database vs. 


Event rate 
/sec (cut) 

Hours of 
Data 

3-event burst 

4-event burst 

exp 

obs 

g 

exp 

obs 

G 

0.0-0. 2 

1759 

62 

73 

+0.95 

0 

2 

1.41 

0.2-0.4 

394 

230 

232 

+0.09 

8 

5 

-0.83 

> 0.4 

64 

236 

257 

+0.95 

12 

5 

-1.70 

Total 

2217 

528 

562 

+1.03 

20 

12 

-1.41 


5 Upper Limit to Density of Exploding PBH 

Based on the null results obtained in this search an upper limit to the local 
density of exploding PBHs can be calculated using the flux of TeV photons 
predicted by the standard model and presented in [12]. 

Previous results of searches for exploding PBH by the Whipple Collabora- 
tion used extrapolated values of collection area and solid angle [18], [16], [23], 
[6], The work presented here involves a search over a larger archival database 
(2217 hours) for bursts on a longer time-scale (1 s) with a lower energy threshold 
(0.4 TeV); hence the sensitivity should be greater. In addition the sensitivity is 
calculated more rigorously than in previous experiments. In the earlier Whip- 
ple experiments the sampling distance was calculated from the trigger efficiency 
rather than gamma-ray collection area. Our new Monte Carlo simulations in- 
dicate that in calculating limits from data taken after 1988 (when the imaging 
system was introduced), the sampling distance was overestimated. In addition, 
a constant gamma-ray sensitivity over the field-of-view was assumed and the 
sensitive volume integrated the estimated sampling distance over the solid an- 
gle of the camera. A more realistic sensitive volume based on Monte Carlo 
simulations is used in this search. Table 3 shows the sampling distances and 
sensitive volumes for a 3- photon burst in the 10m reflector for source offsets out 
to 2° given the gamma-ray flux calculated in [12]. The second column shows 
the collection area for a source at a particular point in the camera and column 
3 shows the maximum distance at which a PBH can be detected. The volume 
probed, V , is calculated by integrating over the solid angle, ftj, covered by the 
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Table 3: Sensitivity of 10m Reflector for 0.4 TeV 7 after applying width and 
leng th cuts at source offsets out to 2°. 


Source Offset 
from Center (°) 

Collection 
Area (10 4 m 2 ) 

Sampling 
Distance (pc) 

Sensitive 

Volume (xl0~ 6 pc 3 ) 

0.00 

5.39 ± 0.89 

0.38 ± 0.15 

1.09 ±0.07 

0.25 

5.07 ± 0.89 

0.37 ± 0.15 

3.03 ±0.20 

0.50 

3.69 ± 0.95 

0.32 ± 0.16 

3.27 ± 0.41 

0.75 

2.73 ± 0.94 

0.27 ± 0.16 

2.75 ± 0.57 

1.00 

2.14 ±0.82 

0.24 ± 0.15 

2.48 ±0.61 

1.25 

1.52 ± 0.78 

0.20 ± 0.14 

1.75 ± 0.60 

1.50 

1.15 ±0.62 

0.18 ± 0.13 

1.51 ± 0.57 

1.75 

0.70 ± 0.48 

0.14 ±0.11 

0.82 ± 0.40 

2.00 

0.31 ± 0.37 

0.09 ± 0.10 

0.25 ± 0.34 

Total sensitive volume | 

16.95 ± 3.77 


annulus (or circle for a source at the centre) described by the source offset: 
y. = This volume is nearly 8 times smaller than the volume 

estimate in Nolan et al. [16] although the hardware system is identical and the 
analysis and burst search methods used in this work are more sensitive. 

The corresponding distances and volumes accessible to the reflector can be 
calculated for the detection of bursts of 4 photons in 1 s. The excess (or deficit) 
of bursts obtained over expectation in the sensitive volume in each of the three 
data groups in the archive is used to find the 99.9% maximum likelihood upper 
limit to the frequency of exploding PBH per year per cubic parsec (Table 4). 
More sensitive upper limits are obtained from air shower experiments which 
operate at higher energies and have longer exposure times; these are shown in 
Table 5. 


Tab! 


e 4: Upper Limits to Exploding PBH Frequency from Archival Search 


Hours 
of Data 

3-event Burst 

4-event Burst 

Volume 
xlO -5 pc 3 

Limit 

xlO 6 pc” 3 yr _1 

Volume 
xlO -5 pc 3 

Limit 

xlO 6 pc~ 3 yr _1 

1759 

394 

64 

1.67 ±0.38 

14.13 ±3.3 
96.2 ± 23.0 
745 ± 170 

1.07 ±0.26 

3.0 ±1.0 
20.4 ± 5.3 
117 ±30 


Models other than those derived from standard elementary particle theory 
have also been invoked to describe the final stages of evaporation of a PBH. 
Those propose a faster process with more degrees of freedom than the standard 
model - the Hagedom model, in which the number of degrees of freedom in- 
creases exponentially with the number of emitted particles, predicts the most 
catastrophic explosion, of 6.0 x 10 34 erg, lasting only 10“ 7 seconds, and com- 
prised mainly of 250 MeV photons. A larger volume is accessible to searches 
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Table 5: Standard elementary particle upper limits to PBH density. 


Experiment 

Energy 

(TeV) 

Reference 

Limit 

pc-V” 1 

CYGNUS 

50 

Alexandreas et al. 1993 [2] 

BBI 

AIROBICC 

20 

Punk et al. 1995 [11] 

ESI 

Tibet 

10 

Amenoruori et al. 1995 [3] 

H 


for this type of event than for PBH exploding according to the standard model. 
Upper limits to Hagedorn-type PBH explosions are, therefore, more restrictive, 
and are given in [17] and [8]. 

6 Conclusions 

A method is described of using an atmospheric Cherenkov telescope to search for 
gamma-ray bursts on short time-scales. The application of this method to the 
search for counterparts to classical gamma-ray bursts detected by the BATSE 
experiment is described elsewhere [7]. Here it is applied to a search through 
a four year database accumulated by the Whipple telescope in its routine dis- 
crete source observing program. The minimum detectable fluence (6 x 10 -9 
erg-cm -2 in one second) is comparable with that of BATSE at much lower en- 
ergies and compares favorably with all other gamma-ray experiments currently 
in operation. 

Null results are obtained and an upper limit to the PBH density is derived, 
where the PBH is exploding via the standard model. This limit of 3.0± 1.0 x 10 6 
pc -3 yr -1 is better than that obtained from previous searches using the Whipple 
10m reflector. It is not, however, as stringent as the PBH limits obtained with 
wide field air shower experiments. The low duty cycle and small field-of-view 
of atmospheric Cherenkov telescopes makes them less efficient than air shower 
experiments like the MILAGRO water-Cherenkov telescope [4] which have a 
large field-of-view and longer exposures. Telescopes with rapid slew speeds are, 
nevertheless, well-suited for the detection of very high energy counterparts to 
BATSE-type bursts [7]. 
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Figure 1: Arrangement of photomultiplier tubes in focus box of 10m reflector. 
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Figure 2: Arrival time difference distribution of successive events in a 60 minute data 
file. 
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Figure 3: 10m image (a) before and (b) after flat-fielding, and (c) after cleaning. 



Figure 4: Representation of Image Parameters. 
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Figure 5: Width (left) and Length (right) distributions for simulated 0.4 TeV photon 
from the zenith (top) and real zenith (bottom) showers. 
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Figure 6: Percentage of events which trigger 10m reflector as a function of impact 

parameter for source offsets (a)0°, (b)0.75°,(c)1.5° for 0.4 TeV showers. 



Figure 7: Trigger efficiency of 10m reflector as a function of source offset and energy. 
An efficiency of 1 is that of an ideal detector which has a 100 % trigger rate for showers 
of impact parameters between 10 and 150 metres, and no triggers outside this range. 
This is equivalent to a collection area of 7 x 10 4 m 2 . 
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Figure 8: Percentage of events which survive Supercuts shape cuts as a function of 
impact parameter for source offsets (a)0 o } (b)0.75 0 ,(c)1.5° for 0.4 TeV showers. 



Figure 9: Ability of 10m reflector to detect gamma rays as a function of source offset 
and energy. Width and Length cuts are made on all events that trigger the system. 
An efficiency of 1 is that of an ideal detector which has a 100% trigger and image 
selection rate for showers with impact parameters between 10 and 150 metres and 
none outside. This is equivalent to a collection area of 7 x 10 4 m 2 . 
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Figure 10: Efficiency curve from simulations (solid line) compared to gamma-ray 
rate measured from Crab Nebula (x) as a function of source offset. 
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Figure 11: Illustration of the most likely points-of-origin (x) along the major axes 

of two events (ellipses). The shaded circles indicate the uncertainty associated with 
these positions. 
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Figure 12: Efficiency of point-of-origin determination from d = 2 — 2 x width /length 
along the major axis. The tolerance on the r-axis is the value by which this calculated 
point-of-origin is allowed to vary from the actual source position. 
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Figure 13: Variation of event rate with elevation raw event rate (top) and event rate 
after selection via image width and length cuts (bottom). 


23 






daqre« (y) 


1.5 

1.0 

0.5 

0.0 

- 0.5 

- 1.0 

- 1.5 


Figure 14: Event centroid distribution of events in a 60 minute data file. 
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Figure 15: Two bursts of four events from real data in 1-second with a possible 
common point-of-origin. 
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Abstract 

Although atmospheric Cherenkov telescopes have restricted fields of view, their fluence sensitivity warrants a search for 
gamma-ray burst phenomena, A search for 400 GeV gamma-ray bursts on a 1 s time-scale using archival data taken between 
1988 and 1992 with the Whipple Collaboration 10 m reflector is presented. No evidence of such bursts is found. Bursts of 
TeV gamma rays have been predicted from exotic astrophysical objects such as Primordial Black Holes and Cosmic Strings. 
An upper limit to the number density of exploding PBH of 3.0± 1 .0 x 10 6 pc“ 3 yr _I is calculated. © 1998 Elsevier Science 
B.V. 

PACS: 95.75.-z; 95.85.Pw; 98.70.Rz; 98.62.Nx 

Keywords: TeV astronomy; Gamma-ray bursts; Primordial black holes 


1. Introduction 

The detection of gamma-ray bursts with peak lu- 
minosity at MeV energies has been one of the most 
exciting discoveries in high-energy astrophysics [9]. 
Although still unexplained after 30 years of intense 
research, the phenomenon has inspired searches for 
counterparts at a variety of wavelengths [14]. It has 


* Corresponding author. ES84, Marshall SpaceFlight Center, AL 
35812, USA; E-mail: vc@msfc.nasa.gov. 


also opened the possibility that there might be similar 
phenomena that have peak luminosity at other wave- 
lengths; here we consider a TeV search for such phe- 
nomena on time-scales of one second using the Imag- 
ing Atmospheric Cherenkov Technique. 

The Whipple Observatory Gamma Ray Collabo- 
ration has been involved in searches for gamma-ray 
bursts since 1978; initially, multi -element first genera- 
tion systems over long baselines were used [ 17] . The 
development of the Imaging Atmospheric Cherenkov 
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Fig. 1. Arrangement of photomultiplier tubes in focus box of 10 m 
reflector 


Technique has significantly improved the sensitivity 
for the detection of point sources [25]. The, Whip- 
ple 10 m imaging camera has successfully detected 
TeV emission from several point sources - the Crab 
Nebula [24], Markarian 421 [19], and Markarian 
501 [20] . It has been the most sensitive instrument of 
its kind for source-centered observations but its off- 
axis sensitivity has not so far been fully character- 
ized [1]. 

Here the response of the Whipple 10 m reflector 
to non-source-centered gamma rays is assessed us- 
ing shower simulations, and a method is developed to 
search for gamma-ray bursts of one second duration. 
We report on a search through 4 years of the Whipple 
database. 

There are a number of exotic suggestions that justify 
a search for gamma-ray bursts in this parameter space. 
These include emission from the decay of primordial 
black holes (PBHs) [13] and cosmic strings [22]. 
The cosmological and physical importance of PBHs 
is well established [12]. As the PBH evaporates, its 
temperature increases in proportion to the mass loss. 
Most models predict a final explosion of energy when 
all possible evaporation channels are available, but the 
number of degrees of freedom of emission is a highly 
controversial issue, and ranges from the conservative 
standard elementary particle model, where it reaches 
a maximum once the three generations of fermions 



Fig. 2. Arrival time difference distribution of successive events 
in a 60 minute data file (histogram). The solid curve shows the 
expected Poisson distribution for this average event rate and total 
number of events. 

are free to be produced, to the runaway Hagedom 
model where the number of emission modes increases 
exponentially beyond this [ 12] . For the conservative 
model, an optimum time-scale of 1 s is found for ob- 
serving the explosion of a PBH using the Whipple 
10 m telescope. During this final second 9.1 x 10 28 
photons are emitted above 0.4 TeV. 

The fluence sensitivity of the detector over this time- 
scale is of order 10 -9 erg cm -2 which is compara- 
ble to or better than that achieved in most bands of 
the gamma-ray spectrum. The sensitivity of the tech- 
nique is compared to that achieved in searches by other 
gamma-ray experiments, operating at energy thresh- 
olds higher than the Whipple instrument. 


2. The Whipple 10 m telescope 

Situated at an altitude of 2.3 km on Mt.Hopkins in 
Arizona, the telescope operated by the Whipple Col- 
laboration has been used as an imaging device since 
1982; the camera has been improved from a 37 to 
a 109 element imaging system during that time. The 
telescope consists of a 10- metre dish (with 248 front- 
coated mirrors) on an alt-azimuth mount and a cam- 
era containing 109 photomultipliers (PMTs). The in- 
ner 91 PMTs are 1.1 cm diameter tubes (Hamamatsu 
R1398) and are connected via amplifiers to trigger 
discriminators. An outer ring of 1 8 tubes is indepen- 
dent of the trigger logic and contributes only to the 
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Fig. 3. 10 m image (a) before and (b) after flat-fielding, and (c) after cleaning. 


imaging of the event. Until 1993, these outer tubes 
were 5 cm tubes giving a total field-of-view of 3.75°. 
They have been replaced by the smaller tubes so that 
the camera field-of-view is now 3°. The arrangement 
of PMTs in the focus box prior to 1993 is shown in 
Fig. 1 and a detailed description of the instrument can 
be found in [5]. 

The camera is triggered if the signal in two or more 
PMTs exceeds a preset threshold in a 10 nsec time 
interval. When the system is triggered, i.e. an event 
is registered, the light level in each PMT is recorded 
over a 25 nsec gate width. The trigger rate during the 
epoch of the observations reported here was typically 
8 events per second. The time is derived from a GPS 
clock and is recorded with each event with an accuracy 
of 0.1 msec. The deadtime is 1.5 msec and the distri- 
bution’ of time differences of successive events from 
a 60 minute data file taken at a telescope elevation of 
80° is shown in Fig. 2 (histogram). It is overlaid on 
the expected Poisson distribution for the same mean 
event rate (solid curve).’ 

The digitized images of the Cherenkov light are ana- 
lyzed off-line. A typical image recorded by the present 


Length 



Fig. 4. Representation of image parameters. 


camera is shown at various stages of processing in 
Fig. 3. Moment-fitting routines are used to define each 
image by an ellipse whose parameters define the event 
(shown in Fig. 4). The event is accepted as a poten- 
tial gamma-ray event or rejected as a hadron-initiated 
event depending on the values of these parameters. 
The axis of the image is defined by minimizing the 
signal-weighted sum of the squares of the perpendic- 
ular angular distances of the pixels. The r.m.s. spread 
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of light perpendicular to and along this axis defines 
the semi-minor and semi-major axes of the ellipse and 
these are known as the width and length of the im- 
age. Image compactness is defined by concentration , 
which is the fraction of the sum over all nonzero pix- 
els of the cleaned image in digital counts (known as 
size) contained in the two or three brightest pixels. 
The other parameters shown in Fig. 4 relate to the ori- 
entation of the image axis relative to the source posi- 
tion in the field-of-view. 

If the position of the source is known (usually at the 
centre of the field of view), then the orientation of the 
image relative to that direction is important for source 
detection. A combination of shape and orientation pa- 
rameters has been used by the Whipple Collaboration 
to reject 99.7% of recorded background while keeping 
50% of gamma rays; the set of selection parameters 
that are most useful for observations of a point-source 
in the centre of the field of view are called Supercuts 
and are defined elsewhere [21 ] . A comparison of the 
width and length domains for photons (simulated) 
and real background (data taken with the telescope 
pointed at the zenith) can be seen in Fig. 5. 

The success of Supercuts and the Imaging Atmo- 
spheric Cherenkov Technique in the detection of point 
sources owes much to Monte Carlo simulations of pho- 
ton and background hadron-initiated cascades. Ow- 
ing to the impossibility of testing the atmospheric 
Cherenkov detector at an accelerator, analytical mod- 
els and simulation of phenomena are used to optimize 
instrumental design and in the development of analy- 
sis techniques. Data taken during observations of the 
Crab Nebula, which plays the role of a TeV gamma-ray 
standard candle, are used to perfect new techniques. 

3. Off-axis sensitivity of 10 m reflector 

The characteristics of the 10 m reflector in point- 
source mode of operation are well understood and de- 
fined [5]: What is less well characterized is the sen- 
sitivity of the instrument when an area of the sky is 
probed for a source without knowing its exact posi- 
tion, i.e., where each point in the field-of-view is a 
potential source. It has been, demonstrated experimen- 
tally that the system is sensitive to a point source of 
gamma rays up to 1° off-axis [ 1,10] . In the work re- 
ported here, the response over the full field-of-view of 


the 10 m reflector is investigated, using simulations, 
in order to use the instrument as a relatively wide field 
burst detector. 

5.7. The off-axis triggering and imaging of gamma 
rays 

The geometrical field-of-view of the 10 m camera 
between 1988 and 1992 was 3.75°, but the imaging 
technique allows reconstruction of images whose ori- 
gins lie outside this area. Conversely, some of the im- 
ages of showers from sources on the edge of the cam- 
era face may fall outside the field-of-view so that the 
camera has reduced collection area for such sources. 

In order to assess the sensitivity of the instrument 
for sources lying at various distances from the centre 
of the camera, a database of nearly 8000 simulated 
gamma-ray-induced showers of energies ranging from 
0.2 TeV to 0.8 TeV was compiled using a simulation 
program based on the KASCADE code of Kertzmann 
and Sembroski [15]. The Whipple collaboration has 
used several different simulation codes in the various 
research institutions (e.g. MOCCA at the University 
of Leeds) and comparison of the simulations shows 
that they are in agreement with each other to within 
10% in the Cherenkov light yields and image param- 
eter distributions that they produce. 

The response of the instrument to the simulated 
gamma-ray showers is gauged by mapping each 
Cherenkov photon in a shower from the mirrors onto 
the imaging camera and building up an event for each 
shower. 

Displacing the focus box from its position at the 
focus of the reflector and constructing events as they 
appear to this displaced camera is equivalent to sim- 
ulating the response of the camera to sources at an 
offset from the centre equal to the displacement of the 
camera but in the opposite direction. The response of 
the camera to gamma-ray showers from sources with 
0° to 2.5° offsets from the centre was investigated by 
successively displacing the camera at 0.25° intervals 
along an arbitrary line, 17° from an axis through the 
centre of a row of PMTs. A different source displace- 
ment angle may be chosen without affecting the re- 
sults. The treatment of simulated images is identical 
to that subsequently applied to real data. The parame- 
terization and cutting procedures execute rapidly rel- 
ative to the shower and instrument simulations, and it 
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Fig. 5. Width (left) and length (right) distributions for simulated 0.4 TeV photon from the zenith (top) and real zenith (bottom) showers. 


was possible to experiment with a variety of parame- 
ters and analysis methods before deciding on a suit- 
able gamma-ray burst search method. 

As the source position moves further away from 
the centre of the camera, less Cherenkov light falls on 
the detector at a given impact parameter. The software 
trigger is a requirement that two of the inner 91 tubes 
register at least 40 photo-electrons and the percentage 
of events triggering the system at each impact param- 
eter is shown in Fig. 6 for source offsets of (a) 0°, 
(b) 0.75°, and (c) 1.5°. 

By weighting the fraction according to the impact 
parameter value (the collection area is bigger at larger 
impact radii), a value for the trigger efficiency over 
the field of view, equivalent to the effective area of the 
detector for each source offset, is calculated, 

.j. . ^ . £, 2vri4>i Sr 

Trigger efficiency = — 1 - — ^ — , ( 1 ) 

where r* ( = i ) is the radius of the annulus at impact 
parameter /, thickness Sr = 2 metres, and <f>i is the 
fraction of the 20 showers at radius r, which would 
trigger the telescope. The summing integer i is incre- 
mented from 12 to 200 in steps of 2. At impact pa- 
rameters below 10 metres, the showers tend to saturate 


the central tube, and the Cherenkov light pool falls off 
rapidly beyond 150 metres, so that we define an ideal 
detector as one which is triggered by all showers with 
impact parameters between 10 and 150 metres, and 
which fails to register any events falling outside that 
area. The trigger efficiency is divided by 7 x 10 4 m 2 
to give the efficiency relative to this arbitrary ideal 
detector. An efficiency of greater than 1 is, therefore, 
possible since showers lying beyond 150 metres may 
trigger the system. This procedure was implemented 
for all 4 shower energies, and the results are shown in 
Fig. 7. Above threshold ( > 0.4 TeV) the instrument 
would register over half of the gamma-ray showers 
from sources out to 1° away from the centre, but it is 
obviously more efficient at detecting those nearer the 
centre than at large offsets. 

Image parameters have traditionally been used, in 
some combination of shape and orientation cuts, to 
discriminate against the hadronic background. In this 
search, however, only the shape parameters are rele- 
vant since there is no preferred direction due to the un- 
known position of the source in the field-of-view. The 
most successful shape selection was found to involve 
cuts in length and width , using the values derived for 
Supercuts, 
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Imped (metres) 


Fig. 6. Percentage of events which trigger the 10 m reflector as a function of impact parameter for source offsets (a) 0°, (b) 0.75°, (c) 
1 .5° for 0.4 TeV showers. 



Fig. 7. Trigger efficiency of the 10 m reflector as a function of 
source offset and energy. An efficiency of 1 is that of an ideal 
detector which has a 100% trigger rate for showers of impact 
parameters between 10 and 150 metres, and no triggers outside 
this range. This is equivalent to a collection area of 7 x 10 4 m 2 . 

0.073° < width <0.15°, 

0.16° < length < 0.30° . 

Fig. 8 shows the percentage of 0.4 TeV events at off- 
sets (a) 0°, (b) 0.75° and (c) 1.5° which trigger the 
camera and satisfy these image shape requirements. 
Combining the percentages for each impact parameter 
produces the results which are shown in Fig. 9. Above 
the detector energy threshold, the collection area for 


gamma rays is fairly flat, and showers from sources 
which are up to half a degree off-axis are still selected 
over 50% of the time. There are two competing ef- 
fects which explain the shape of the curves: at high 
energies, showers with large impact parameters pro- 
duce enough Cherenkov light to trigger the telescope 
more often than distant low-energy showers, but the 
Supercuts selection is biased towards small showers, 
so that the higher energy showers are less likely than 
the low-energy events to survive the width and length 
cuts regardless of impact parameter. At large source 
offsets the fraction of events selected decreases with 
shower energy even though the high-energy events ap- 
pear to trigger the camera more often than at lower 
energies. If the showers triggering the system at these 
large offsets are examined, it can be seen that they are 
the longest images in the high-energy databases. 

Observations of the Crab Nebula at angular offsets 
from the centre of the field-of-view of the 10 m tele- 
scope were used to confirm the efficiency curves de- 
rived from these simulations. The Crab Nebula is a 
steady source of TeV gamma rays, and by comparing 
the rate of gamma rays measured after applying the 
Supercuts selection technique to data from September 
1996 at offsets ranging from 0° to 1.5°, one obtains 
the efficiency for gamma ray collection displayed in 
Fig. 10. The errors in Fig. 10 reflect the variation in 
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Fig. 8. Percentage of events which survive Supercuts shape cuts as a function of impact parameter for source offsets (a) 0°, (b) 0.75°, 
(c) 1.5° for 0.4 TeV showers. 



Fig. 9. Ability of 10 m reflector to detect gamma rays as a function 
of source offset and energy. Width and Length cuts are made on all 
events that trigger the system. An efficiency of 1 is that of an ideal 
detector which has a 100% trigger and image selection rate for 
showers with impact parameters between 10 and 150 metres and 
none outside. This is equivalent to a collection area of 7 x 10 4 m 2 . 

evenf rates for the offset observations, and are large 
owing to the small exposure time: only 1 hour of 
source observations was made at a source offset of 
0.5° from the centre of the field-of-view, 1.5 hours at 
1.0°, and 2 hours at 1.5° offset. Because the field-of- 
view of the telescope was smaller when these offset 
observations were made compared to that simulated 
in this study, a perfect correlation is not expected, and 



Fig. 10. Efficiency curve from simulations (solid line) compared 
to the gamma-ray rate measured from Crab Nebula (r) as a 
function of source offset. 

one anticipates that the efficiency at larger source off- 
sets be slightly lower relative to source-centered ob- 
servations than in the simulations. 

5.2. Orientation of gamma- ray images 

The orientation of the ellipse fitted to each image 
is represented by its major axis, and the most likely 
point-of-origin of the shower progenitor on the field- 
of-view lies on this axis at a distance d in degrees 
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Fig. 11. Illustration of the most likely points-of-origin ( x ) along 
the major axes of two events (ellipses). The shaded circles indicate 
the uncertainty associated with these positions. 

related to the ellipticity of the image, 

d = 2 — 2(width I length) . ' (2) 

This algorithm was developed based on simulations 
of gamma-ray showers from sources at the centre of 
the field-of-view, and was found to be accurate to about 
0.3° either side of this point [ 1 ] for source-centered 
observations. The relation between the image, its cen- 
tre, and the most likely point-of-origin are shown in 
Fig. 11. 

Because the ellipses derived from the moment- 
fitting routines are symmetrical, the point-of-origin 
may lie either side of the centre of the image. The 
following equations are used: 

\/ ( y’or >’cen ) ^ "b ( ^or -^cen ) ^ 

' " = 2 — 2 x width/length , 
y or = m x x 0T + c , 

to find the most likely points-of-origin x OI , y or on both 
sides of the centre of light .t ccn . y C en along the major 
axis y = m x x + c. The distance 

\/(y0 - >’or) 2 + Uo - *or) 2 , (3) 

where xo. yo are the true positions of the simulated 
source, is a measure of how well Eq. (2) finds the 
point-of-origin of a shower. The tolerance is defined 



Fig. 12. Efficiency of point-of-origin determination from 
d = 2 — 2 x width/ length along the major axis. The tolerance on 
the *-axis is the value by which this calculated point-of-origin is 
allowed to vary from the actual source position. 


as the maximum allowed value of this distance, in 
degrees. A range of tolerance levels from 0.1° to 
0.7° was explored to determine the off-axis efficiency 
of the point-of-origin determination, and the results 
are displayed in Fig. 12. The total efficiency & for 
each energy was obtained by combining source po- 
sitions and weighting each efficiency (c fi r ) according 
to source offset (r). The efficiency values featured 
are relative to an ideal detector with <f> r = 1 at each 
source offset where the gamma-ray images are se- 
lected and their origin is correctly determined so that 
<P = <f> r r/ r, with r varying from 0.25° to 2.5°. 
The efficiency curve is fairly constant as the source is 
moved away from the centre and appears independent 
of impact parameter and photon energy. 

Real background data were subjected to the same 
analysis so that the effect of varying the point-of-origin 
tolerance value on background acceptance could be 
assessed. It was found that in looking for bursts of 3 
or more events over a 1 -second time-scale, the highest 
signal to noise ratio could be obtained with a toler- 
ance level of 0.425°. This selection criterion enables 
66 ± 5% of shape-selected simulated events (over the 
impact parameter range 12 to 200 metres and source 
offset 0° to 2.5°) to be correctly located. The shaded 
areas in Fig. 1 1 show the regions from which the im- 
aged events could have originated with this tolerance 
value. 
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Fig. 13. Variation of event rate with elevation raw event rate (left) and event rate after selection via image width and length cuts (right). 


It is then possible to apply this calculated off-axis 
response of the detector to the search for counterparts 
to the delayed component of BATSE bursts (reported 
elsewhere [7] ) and to the serendipitous detection of 
gamma-ray bursts during the normal operation of the 
gamma-ray telescope (described below) ^ 

4. Searches for 1 s TeV gamma-ray bursts 

The archival data of the Whipple Observatory taken 
with the 10 m reflector with an energy threshold above 
0.4 TeV between September 1988 and September 1992 
were used for this search. In general these observa- 
tions were in the form of files of (approximately) 28 
minute duration containing information on each event 
(the arrival time, the outputs of the 109 ADCs from 
the 109 pixels, and other housekeeping data). Exami- 
nation of the observing logs allowed the identification 
and rejection of 1% of the data including (i) obser- 
vations made when there were instrumentation prob- 
lems, (ii) data taken with discriminator thresholds or 
high-voltage values significantly different from those 
used in normal operation, and (iii) observations made 
white testing new or experimental configurations (for 
example, experiments with filters). A further 1% of 
the files were rejected because of obvious anomalies, 
such as a large number of bright stars in the field-of- 
view, or because they, were taken at a telescope ele- 
vation less than 35° above the horizon. Following re- 
jection of unsuitable data, 2217 hours of observations 
comprised the database. The database was divided ac- 


cording to the event rate of candidate gamma rays. 
Variation in event rates was caused by differences in 
zenith angle, mirror reflectivity and discriminator set- 
tings. The elevation of the source under observation is 
known to affect the energy threshold of the detector, 
and also influences the appearance of both gamma and 
hadronic images. A larger depth of atmosphere must 
be penetrated by the extensive air shower particles if 
they are to produce light lower in the atmosphere, so 
that a telescope pointed close to the zenith is sen- 
sitive to lower energy showers than if it is operated 
nearer the horizon. Consequently, the raw event rates 
decrease with increasing zenith angle - this relation- 
ship is shown in Fig. 13 (left) . Despite this overall fall 
in event rates at low elevations, the Cherenkov image 
on the camera face is narrower on the horizon than 
at the zenith, and this effect results in more hadronic 
events being selected as gamma rays on the basis of the 
width parameter. Fig. 13 (right) shows the variation 
of image-selected event rate with telescope elevation. 

Nearly 80% of observations were made with the 
telescope elevation above 55°, where data rates and 
image characteristics are fairly homogeneous. These 
observations form the bulk of the 1759 hours of data 
which had an average rate of < 0.2 candidate gamma- 
ray events per second, and comprise the largest of the 
three subsets of the Whipple database. The rest of the 
data was divided according to the average event rate; 
rate = 0.2 to 0.4 per second (394 hours) and rate > 
0.4 events per second (64 hours). Standard routines 
were used to flat-field the data, tubes with bright stars 
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Rg. 14. Event centroid distribution of events in a 60 minute data 
file. 

in their field-of-view were turned off in hardware or 
software, and the image parameters were calculated 
as described in [21]. The events were then charac- 
terized by these image parameters (and their arrival 
times) . The evenness of the camera response across 
the field-of-view is illustrated in Fig. 14, showing the 
distribution of event centroids in camera coordinates 
(degrees from centre) for a 60 minute observation. 

4.L The burst search 

Initially the events were culled to include only those 
events whose images satisfied the gamma-ray event 
criteria. Events which survived the width and length 
cuts described above were included in a reduced data 
set which was submitted to a search for bursts of 3 or 
more events in a 1 s time interval. 

. .The arrival time of the first event was taken as the 
possible start of a 1 $ burst and was compared to the 
time of the second event. If the time difference was un- 
der 1 s, the next event was examined. If not, the start of 
the time window was moved to the arrival time of the 
second event. The window was moved through the file 
until all events in a file have been processed, and the 
total number of bursts with from 3 to 10 shape- selected 
events in the file is obtained. No candidate bursts were 
found with more than 4 events. The number of 3- and 
4-fold burst events are shown in Table 1 in the col- 
umn headed “obs”. A background (control) file was 


generated from each real data file by scrambling the 
events in each parameterized file, while maintaining 
the original time sequence. When the background file 
was subjected to the burst-search program, the image- 
selected events were the same as those picked from 
the real data file, but in a random order and attached to 
the original arrival times of other events. The number 
of 3- and 4-fold bursts thus found is shown in Table 1 
headed “exp”. The difference between the observed 
and the expectation is expressed in terms of the stan- 
dard deviation, cr defined as (obs -b exp) 1 / 2 . In no 
case is the difference statistically significant. 

To this point no use has been made of the fact that 
the events in the putative 3- and 4-fold bursts must 
come from the same point in the sky. Hence all the 
bursts were subjected to the common origin selection 
described above. Less than 2% of the putative bursts 
passed this selection. These are listed as before in 
Table 2; the expectation is derived by applying the 
common origin selection to the scrambled events. 

Again, no excess of 3- or 4-fold bursts of candidate 
gamma-ray events on a 1 -second time-scale is found 
in the 2217 hours of the Whipple database. Two can- 
didate four-fold events are shown in Fig. 15. If a TeV 
component to a background of bursters exhibiting ac- 
tivity on a 1 s time-scale exists, then the number of 
these bursts is low and below the sensitivity of the in- 
strument in its current form. 


5. Upper limit to density of exploding PBH 

Based on the null results obtained in this search, an 
upper limit to the local density of exploding PBHs can 
be calculated using the flux of TeV photons predicted 
by the standard model and presented in [ 12] . 

Previous results of searches for exploding PBH by 
the Whipple Collaboration used extrapolated values of 
collection area and solid angle [ 18,16,23,6] . The work 
presented here involves a search over a larger archival 
database (2217 hours) for bursts on a longer time- 
scale ( 1 s) with a lower energy threshold (0.4 TeV); 
hence the sensitivity should be greater. In addition the 
sensitivity is calculated more rigorously than in previ- 
ous experiments. In the earlier Whipple experiments 
the sampling distance was calculated from the trigger 
efficiency rather than gamma-ray collection area. Our 
new Monte Carlo simulations indicate that in calculat- 
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Table 1 


Frequency of gamma-ray burst candidates in Whipple database vs expected burst frequency (no common origin sought) 


Event rate/sec (cut) 

Hours of data 

3-event burst 


4-event burst 



exp 

obs 

(T 

exp 

obs 

cr 

0.0-0.2 

1759 

4747 

4762 

+0.15 

234 

245 

+0.50 


0.2-0.4 

394 

13463 

13679 

+ 1.31 

1308 

1286 

-0.43 


> 0.4 

64 

10597 

10539 

-0.40 

1800 

1724 

-1.28 


Total 

2217 

28807 

28980 

+0.72 

3342 

3255 

-1.07 


Table 2 









Frequency of gamma-ray burst candidates in Whipple database vs expected burst frequency (with common origin) 



Event rate /sec (cut) 

Hours of data 

3-event burst 


4-event burst 





exp 

obs 

cr 

exp 

obs 

a 


0.0-0.2 

1759 

62 

73 

+0.95 

0 

2 

1.41 


0.2-0.4 

394 

230 

232 

+0.09 

8 

5 

-0.83 


> 0.4 

64 

236 

257 

+0.95 

12 

5 

-1.70 


Total 

2217 

528 

562 

+ 1.03 

20 

12 

-1.41 







Fig. 15. Two bursts of four events from real data in 1 second with a possible common point-of-origin. 


ing limits from data taken after 1988 (when the imag- 
ing system was introduced), the sampling distance 
was overestimated. In addition, a constant gamma- 
ray sensitivity over the field-of-view was assumed and 
the sensitive volume integrated the estimated sampling 
distance over the solid angle of the camera. A more 
realistic sensitive volume based on Monte Carlo sim- 
ulations is used in this search. Table 3 shows the sam- 
pling distances and sensitive volumes for a 3-photon 


burst in the 10 m reflector for source offsets out to 2° 
given the gamma-ray flux calculated in [12]. The sec- 
ond column shows the collection area for a source at 
a particular point in the camera and column 3 shows 
the maximum distance at which a PBH can be de- 
tected. The volume probed, V iy is calculated by inte- 
grating over the solid angle, covered by the annu- 
lus (or circle for a source at the centre) described by 
the source offset: Vi = This volume 
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Table 3 

Sensitivity of 10 m reflector for 0.4 TeV y after applying width and length cuts at source offsets out to 2 


Source offset from center (°) 

Collection area (10 4 m 2 ) 

Sampling distance (pc) Sensitive volume (xlO 6 pc 3 ) 

0.00 

5.39 ±0.89 

0.38 ±0.15 

1.09 ±0.07 

0.25 

5.07 ±0.89 

0.37 ±0.15 

3.03 ± 0.20 

0.50 

3.69 ± 0.95 

0.32 ±0.16 

3.27 ±0.41 

0.75 

2.73 ± 0.94 

0.27 ±0.16 

2.75 ± 0.57 

1.00 ... - 

2.14 ±0.82 

0.24 ±0.15 

2.48 ±0.61 

1.25 

1.52 ±0.78 

0.20 ±0.14 

1.75 ±0.60 

1.50 

1.15 ±0.62 

0.18 ± 0.13 

1.51 ±0.57 

1.75 

0.70 ±0.48 

0.14 ± 0.1 1 

0.82 ± 0.40 

2.00 

0.31 ±0.37 

0.09 ±0.10 

0.25 ± 0.34 

Total sensitive volume 



16.95 ±3.77 

Table 4 




Upper limits to exploding PBH frequency from archival search 



Hours of data 3 -event burst 


4-event burst 

volume xlO - 

5 pc 3 limit xlO 6 pc' 

-3y r -i volume x 1 0~ 5 pc 3 limit xl0 6 pc^ 3 yr“ l 

1759 1.67 ±0.38 

14.13 ± 3.3 

1.07 ±0.26 3.0 ±1.0 

394 

96.2 ±23.0 


20.4 ± 5.3 

64 

745 ± t70 


117 ±30 


is nearly 8 times smaller than the volume estimate 
in Nolan et al. [16] although the hardware system is 
identical and the analysis and burst search methods 
used in this work are more sensitive. 

The corresponding distances and volumes accessi- 
ble to the reflector can be calculated for the detection 
of bursts of 4 photons in 1 s. The excess (or deficit) of 
bursts obtained over expectation in the sensitive vol- 
ume in each of the three data groups in the archive 
is used to find the 99.9% maximum likelihood upper 
limit to the frequency of exploding PBH per year per 
cubic parsec (Table 4). More sensitive upper limits 
are obtained from air shower experiments which oper- 
ate at higher energies and have longer exposure times; 
these are shown in Table 5. 

Models other than those derived from standard ele- 
mentary particle theory have also been invoked to de- 
scribe the final stages of evaporation of a PBH. Those 
propose a faster process with more degrees of free- 
dom than the standard model - the Hagedom model, 
in which the number of degrees of freedom increases 
exponentially with the number of emitted particles, 
predicts the most catastrophic explosion, of 6.0 x 10 34 


Table 5 

Standard elementary particle upper limits to PBH density 


Experiment 

Energy 

(TeV) 

Ref. 

Limit 

(pc - 3 yr -1 ) 

CYGNUS 

50 

[2] 

6.1 x 10 5 

AIROBICC 

20 

[11] 

8.9 x 10 5 

Tibet 

10 

[3] 

4.6 x 10 5 


erg, lasting only 10 -7 seconds, and comprised mainly 
of 250 MeV photons. A larger volume is accessible to 
searches for this type of event than for PBH explod- 
ing according to the standard model. Upper limits to 
Hagedom-type PBH explosions are, therefore, more 
restrictive, and are given in [17,8]. 

6. Conclusions 

A method is described of using an atmospheric 
Cherenkov telescope to search for gamma-ray bursts 
on short time-scales. The application of this method 
to the search for counterparts to classical gamma-ray 
bursts detected by the BATSE experiment is described 
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elsewhere [ 7] . Here it is applied to a search through a 
four year database accumulated by the Whipple tele- 
scope in its routine discrete source observing program. 
The minimum detectable fluence (6 x 10“ 9 erg cm -2 
in one second) is comparable with that of BATSE at 
much lower energies and compares favorably with all 
other gamma-ray experiments currently in operation. 

Null results are obtained and an upper limit to the 
PBH density is derived, where the PBH is explod- 
ing via the standard model. This limit of 3.0 ± 1.0 x 
10 6 pc“ 3 yr“ 1 is better than that obtained from pre- 
vious searches using the Whipple 10 m reflector. It 
is not, however, as stringent as the PBH limits ob- 
tained with wide field air shower experiments. The 
low duty cycle and small field-of-view of atmospheric 
Cherenkov telescopes makes them less efficient than 
air shower experiments like the MILAGRO water- 
Cherenkov telescope [4] which have a large field- 
of-view and longer exposures. Telescopes with rapid 
slew speeds are, nevertheless, well suited for the de- 
tection of very high-energy counterparts to BATSE- 
type bursts [7] . 
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ABSTRACT 

Based on BACODINE network notification the Whipple Observatory gamma-ray telescope has been 
used to search for the delayed TeV counterpart to BATSE-detected gamma-ray bursts. In the fast slew 
mode, any point in the sky can be reached within two minutes of the burst notification. The search 
strategy, necessary because of the uncertainty in burst position and limited FOV of the camera, is 
described. 

INTRODUCTION 

The Whipple collaboration has an active observational program which is dedicated to searching for 
TeV counterparts of classical gamma ray bursts. Since May 1994, the Whipple Observatory 10 me- 
ter gamma-ray telescope (Cawley, et al., 1990) has made 16 follow-up observations of BATSE burst 
notifications, but to date no positive identifications have been made (Connaughton, et al., 1995, Con- 
naughton, et al., 1997a). The field of view of our current instrument is limited to 3.5 degrees and the 
error box for the BATSE position of gamma-ray bursts is 10 degrees in diameter. Therefore, more 
than a single pointing has to be used for the TeV follow up after a gamma-ray burst. The current 
observation technique is to take a 28 minute exposure at the initial location and 4 more exposures 3 
degrees away, with a final exposure taken again at the initial position as in Figure 1 a (Connaughton, et 
al., 1997a). In order to find any associated TeV component a two dimensional (2D) search strategy has 
been developed in order to locate the precise region of emission (Akerlof, et al., 1993, Connaughton, 
et al., 1997a, Lessard, et al., 1997). 

Assuming gamma ray bursts have an associated TeV component, it is uncertain as to just how long 
after the BATSE notification of a gamma-ray burst that this emission would last. At present, we cover 
less than 50% of the actual BATSE error box in 2.5 hours and of the 16 follow-ups observed to date, 
only 6 had full observations of all positions (Table 1). All other follow-ups were terminated after 
an hour due to sun/moonrise or the position being too low to track. Owing to the uncertainty in the 
timescale of bursts and the lack of coverage of the whole error box, a non-detection cannot be used 
as an argument against the possibility of gamma-ray burst emission extending up to TeV-emission 
energies. 



Table 1: Whipple Observations of BACODINE Positions 


Date 

BATSE 

Trig. 0 

BACODINE 

Intensity 

Cts/s b 

Difference 
BACO- 
Hunts. (deg) c 

Durat. 

BATSE 

(s) d 

Whipple 

delay 

(min)* 

Elev- 

ation 

(deg/ 

Cycles 

observed 

0-6/ 

940516 

2980 

630 

N/A 

N/A 

19 

24-14 

1 (1 hour) 

950208 

3408 

778 

1.64 

N/A 

16 

N/A 

1,2 

950405 

3494 

704 

9.5 

N/A 

8 

27-33 

1 ,2,3+cont 

950524 

3598 

8726 

1.37 

6 

5 

56-31 

complete 

950625 

3649 

1661 

3.84 

40 

18 

28-41 

1-6 

950701 

3658 

9134 

2.3 

15 

56 

41-69 

complete 

951117 

3909 

1955 

10.1 

25 

5 

30-24 

1 (1 hour) 

951119 

3911 

801 

6.98 

60 

20 

45-54 

complete 

951124 

3918 

1231 

5.98 

150 

2 

59-76 

complete 

951220 

4048 

6698 

2.88 

17 

27 

71-37 

complete 

960521 

5467 

2919 

9.13 


14 

80-68 

1,2 

961017 

5634 

3193 

2.74 

41 

7 

28-23 

complete 

961 111* 

5665 

496 

- 

- 

4 

54-80 

1-4 

961206 

5706 

1312 

6.4 

1.4 

5 

26-22 

1-3 

970304 

6113 

8842 

3.1 

11 

3 

28-42 

1-4 

970308 

6117 

1120 

13.6 

2 

4 

43-72 

complete 


"Trigger number of BATSE burst 

^Number of counts/s measured by BATSE in the first 1 or 2 seconds of the burst. 

c Angular difference between the burst coordinates provided through BACODINE and the final estimated 
burst position from Huntsville. 

‘'Duration of the burst seen by BATSE. 

'Length of time that elapsed between sending the BACODINE message and the start of the first Whipple 
run. 

/ Average elevations of the first and last positions covered by the 10-metre telescope. 

* Whipple coverage of burst area. The numbers listed indicate the positions covered in the observations. 
The 6 positions are described in the text and are shown in Figure la. An entry of ’’complete” means all 6 
positions and their control observations were completed. 

'’Trigger #5665 was a terrestial trigger 


RASTER SEARCH TECHNIQUE 

In order to achieve a better coverage of the 10 degree error box we now scan a circular region around 
the BATSE burst position with a constant angular velocity. A raster scan across parallel lines in 
declination (Figure lb) with a spacing of 1 degree results in complete coverage of the BATSE error 
box. A continuous motor speed of 0. 1 degree/s facilitates a scan of the whole error box within 400 
seconds. The fast scanning method has the advantage that we may give a flux upper limit within 
the first 300-400 seconds or less depending on the final improved BATSE burst position. Although 
the sensitivity within the first 300 seconds is not great, a burst of the calibre of 910503 or 930131 
would be easily detected if the spectrum extends from the GeV region up to TeV energies ( Hurley 
1996). The time scale of those two bursts was 84 and 100 seconds respectively and only if the burst 
positions are close to the actual telescope position can a resonable coverage of the BATSE error box 
be achieved. 

Recent improvements have lowered the time for response to a BACODINE notification of the Whipple 
telescope to within 15 seconds of the BATSE trigger. A new tracking system for the telescope will be 
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Fig. 1: (a): Current Response to BACODINE notifications. Each circle represents the field of view of 
the 10-meter reflector at the 5 different positions following a BATSE burst, (b): Raster Scan response. 
The reflector is moved at a constant velocity continuously along the track to cover the BATSE error 
box. 


implemented in the fall of 1997, making the most distant point on the sky accessible within 2 minutes. 
A burst seen by BATSE and EGRET on 940217 lasted for some 1.5 hours at GeV energies. This 
class of gamma-ray burst might easily have been missed by the previous method but it could now 
be observed with a reasonable sensitivity (9 minutes effective observation time assuming a 3 degree 
FOV and a 10 degree error box). 

The result of each raster scan will be analysed and displayed by a quicklook analysis to search for a 
strong excess. If there is no significant excess this scanning process can be continued for 2-5 hours if 
possible. If something significant is apparent, the telescope will be pointed at this position. 


ANALYSIS 


The Whipple high energy gamma-ray telescope im- 
ages Cerenkov light from air showers onto an array 
of 151 photomultipliers covering a 3.5 degree field 
of view. The imaging technique uses information of 
the angular spread and orientation to reject more than 
99.7% of the cosmic-ray induced showers while re- 
taining over 50% of the possible gamma-ray induced 
events. For discrete objects the source is placed at the 
center of the field of view and candidate gamma-ray 
events (CGRE) are selected on the basis of the a com- 
bination of (a) image shape (length & width ) and (b) 
orientation of the major axis of the shower (Reynolds, 
et al., 1993). However, as a burst source has an un- 
known position a 2-D search can determine the inci- 
dent direction of the CGRE. The search makes use of 
the orienation, elongation and asymmetry of the im- 
age. Monte Carlo studies have shown that gamma-ray 



Fig. 2: 2-D reconstruction of points of ori- 
gin of showers. For small impact parame- 
ters, the image has a structure close to that 
of a circle, whereas with increasing impact 
parameter it appears more elongated. 





images are a) aligned towards the source position b) elongated in proportion to their impact parameter 
on the ground and c) asymetrical, images are weighted towards their point of origin (Figure 2). The 
distance between the image centroid and the true point of origin of a shower is represented by the 
imaging parameter disp: 


disp = 1 .85 - 1 .85 {width/ length ) ( 1 ) 

from which a unique point of origin can be calculated. As each event represents in itself a unique 
telescope position in Right Ascension (RA) and Declination (Dec), each point of origin must be 
translated into an absolute co-ordinate grid. Each point on the grid is then tested for an excess of 
gamma- ray-like events. The off center properties of the Whipple gamma-ray telescope are well es- 
tablished (Connaughton, et al., 1997b) and are used to determine a flux at any point in the camera. 
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ABSTRACT 

The Crab Nebula has become established as the standard candle for TeV gamma-ray astronomy using 
the atmospheric Cherenkov technique. No evidence for variability has been seen. The spectrum 
of gamma rays from the Crab Nebula has been measured in the energy range 500 GeV to 8 TeV 
at the Whipple Observatory by the atmospheric Cherenkov imaging technique. Two methods of 
analysis involving independent Monte Carlo simulations and two databases of observations (1988-89 
and 1995-96) were used and gave close agreement. Using the complete spectrum of the Crab Nebula, 
the spectrum of relativistic electrons is deduced and the spectrum of the resulting inverse Compton 
gamma-ray emission is in good agreement with the measured spectrum if the ambient magnetic field 
is about 25-30 nT. 

1. INTRODUCTION 

The Crab Nebula has become a standard candle in TeV astronomy; it has been detected by many 
groups and its integral flux appears constant. (We refer to the “steady” emission, not emission modu- 
lated at the pulsar period. We have not detected the latter at TeV energies. See paper by G. Gillanders 
et al., in these proceedings.) It is also well on the way to becoming a standard candle with regard to 
TeV spectral content. Synchrotron emission from the Crab covers a remarkably broad range termi- 
nating at about 10 8 eV, where a new component attributed to inverse Compton scattering begins. It is 
this component that we detect. The TeV spectrum is sensitive to the primary electron spectrum, the 
nebular magnetic field and the spatial distribution of electrons and magnetic field within the nebula. 

In this paper we briefly describe two methods for extracting TeV spectra, compare results from 
the Whipple Observatory Imaging Cherenkov Telescope for the 1988/89 and 1995/96 observing sea- 
sons and comment on implications for the physics of the nebula. The methods, developed at Iowa 
State University and the University of Leeds, are based on Monte Carlo simulations using completely 
independent code and use different approaches in determining the overall gain of the Cherenkov tele- 
scope. In the ISU approach, the gains of the photomultiplier tubes, mirror reflectivities, etc., are 
measured and combined to find the overall gain. In the Leeds approach, the observed brightness 
distribution of cosmic-ray images is combined with cosmic-ray simulation results to determine the 
overall gain of the telescope. The methods are described in detail in “Paper I,” Mohanty et al., (1997), 
and the resulting TeV spectrum is put into context of other observations with implications for the 
physics of the nebula in “Paper II,” Hillas et al., (1997). The latter paper also compares our Crab TeV 
spectrum with those of other groups. 



2. METHOD 1 : A TRADITIONAL APPROACH 

A straightforward approach to the determination of TeV spectra was developed at Iowa State Uni- 
versity. Three components are required. First, a method of distinguishing gamma-ray images from 
background cosmic-ray images. The standard method is to use “supercuts" as described in, e.g., 
Punch et al. (1996). The images are characterized by second order moments giving the width, length, 
distance of the image centroid from the optic axis and alpha, the angle by which the image major axis 
misses passing through the optic axis. More that 99% of the background can be rejected by requiring 
small values of width, length and alpha. However, this procedure results in a strong bias against the 
images of higher energy gamma rays which tend to be longer and broader and hence more cosmic-ray 
like. Mohanty (1995) has modified the procedure so that the image selection criteria depend on the 
total brightness or size of the image as well. (The size can be used as an estimate of the energy.) 
The telescope collection areas for simulated gamma rays for the 1995/96 season using standard and 
“extended” supercuts is shown in Fig. 1 . 

The second component needed is 
a way to estimate the energy of each 
gamma-ray image. Two desirable crite- 
ria are (a) good resolution and (b) negli- 
gible bias. The former is important to de- 
tect small structures in the spectrum and 
the latter is important to avoid distortions. 

We obtained a resolution of AE/E ~ 0.36 
with negligible bias by using a second 
order polynomial in size and distance as 
described in Paper I. The energy resolu- 
tion function is, to a good approximation, 

Gaussian in the variable log(£). It is plot- 
ted at an arbitrary energy in Fig. 1 . 

In data taken for spectral analysis, 
each on-source observation is followed by 
an “off-source” observation covering the 
same range of range of azimuth and ele- 
vation angles. The images for both types 
of observation are selected for gamma-like events. The estimated energies from corresponding obser- 
vations are histogrammed and the difference ascribed to gamma rays from the source. This can then 
be fit by a power law or the fluxes extracted as described in Paper I. 

3. METHOD 2: AN INTUITIVE APPROACH 

A different approach with emphasis on verifiability was developed at the University of Leeds. There 
are two pieces to this approach, a method for selecting images likely to have been initiated by cosmic 
gamma rays and a method for determining the primary gamma-ray energy spectrum from the observed 
size spectra. Earlier descriptions of this approach are in (Hillas and West, 1991) and (West, 1994). 

The selection criterion is a “cluster” or “spherical” method in which a single parameter is used 
to characterize the gamma-ray-like nature of an image and correlations between image parameters are 
incorporated naturally. Simulated gamma rays produce images with four parameters (width, length, 
distance and alpha) that populate a four-dimensional space. Each real image can be tested to see if 
it is likely to be a gamma-ray image by the value of the Mahalonobis distance between it and the 
centroid of the cluster. (This is equivalent to scaling and rotation of the axes so that the window is 
spherical and fluctuations uncorrelated.) The selection window is defined by the expected position 
and dimensions of the gamma-ray parameter cluster for several broad ranges of image size. 



1 10 
Log(Energy) 

Fig. 1: The collection areas for standard and extended 
supercuts for the 1995/96 observing season are shown 
above. The resolution function is also plotted. 




In order to extract a spectrum, a size histogram is then computed for on-source and off-source 
observations and the difference histogram is ascribed to gamma rays. A simulated size spectrum is 
then computed starting with a power-law primary energy spectrum for gamma rays. A weight is given 
to each simulated gamma ray and, by adjusting these weights, the spectrum can be varied so that its 
size spectrum matches that of the difference histogram. This method is simple, easy to implement and 
avoids the complexity of calculating collection areas and bias-free energy estimates. It is described in 
detail in Paper I. 

4. RESULTS 


The spectra obtained from the 1988/89 
and 1995/96 seasons using Method 1 and 
1988/89 season using Method 2 are in 
good agreement as shown in Fig. 2. As de- 
scribed in Paper I, we have tested the sen- 
sitivity of the results to uncertainties in the 
Monte Carlo simulations and have found 
that the results are relatively robust. The 
combined TeV fluxes from both seasons 
from Method 1 are well fit with the simple 
power law spectrum in which the differen- 
tial flux (7(E)) is given by: 

(3.3 ± .2 ± .7) io" 7 (-^)- 2 - 45± 08± 05 

(!) 

in units of m -2 s -1 TeV -1 where the first 
errors are statistical and the second are our 
estimate of systematic errors. 

However, the simple extrapolation of 
this fit to lower energies yields fluxes far in 
excess of those observed by EGRET as is 
clear from Fig. 3 which shows a quadratic 
fit of log(7) vs. log(E) to our data and to an 
averaged point representing the EGRET flux (Nolan etal., 1993) at 2 GeV. This fit may be written: 

7(E) = (3.25) 10 _7 (£/TeV) _2 ' 44_a,35log ' 0 ^m _2 s“ 1 TeV _1 . (2) 

5. COMMENTS ON INTERPRETATION 

Most of the early inverse Compton models for TeV gamma rays assumed a constant magnetic field in 
the principal source region where these are produced (e.g., Gould 1965 or Rieke and Weekes 1969) 
whereas more recent models (De Jager and Harding, 1992 or Aharonian and Atoyan) incorporate 
hydrodynamic plasma/field flow making the calculations more complex and the results probably more 
realistic. Here, we try to stay close to the data and make the simpler assumption that the field is 
constant. The broad synchrotron emission band apparently extends up to 10 8 eV and is boosted 
to higher energies via inverse Compton scattering. The scattering giving rise to TeV gamma rays 
occurs in the Klein-Nishina rather than in the Thomson scattering regime. This implies that the 
electrons giving rise to our detected gamma rays must have energies in the range of 2-10 TeV and the 
corresponding scattered photons would mostly have energies of about 0.005 to 0.3 eV. This conclusion 
is only very weakly model dependent (see Paper II). Since electrons with energies of a few TeV 
generate synchrotron radiation at about 0.4 keV in a field of about 25 nT (see next paragraph), the 



Fig. 2: The Crab spectrum in the range 0.3 to about 
8 TeV extracted using Methods 1 ( open circles) and 2 
(solid circles) for the Whipple 1988/89 database and us- 
ing Method 1 (x's) for the 1995/96 database are shown 
above. Also shown is a fit to the combined Method 1 re- 
sults (solid line) as well as an earlier spectrum (dashed 
line) taken from Vacanti etal. (1991). 




Einstein Observatory X-ray pictures of Hamden and Seward (1984) also show the part of the nebula 
emitting TeV gamma rays. 

For assumed magnetic field values 
and the observed synchrotron flux, it is 
possible to deduce the spectrum of pri- 
mary electrons, presumably generated in 
the shock at the termination of the pul- 
sar wind (see, e.g., Coroniti and Kennel 
1985). From the ambient photon density 
and deduced electron spectrum, the TeV 
flux can be calculated, and results for B 
fields of 10, 20 and 40 nT are shown in 
Fig. 3. The shaded regions reflect un- 
certainties arising from ill-defined UV to 
soft-X-ray region of the synchrotron pho- 
ton spectrum. As can be seen from the fig- 
ure, the effective B field must lie between 
20 and 40 nT with 27 nT falling very near mipple and EGRET observations of the Crab 

to our TeV data, mce even more ener g amma . ra y spectrum. The dotted line is a fit to the 
getic electrons only keep their energy for a wm k and a si k 2 GeV flux value. Full-line 

short time, they should exist only near the .... ~ „ a r , . 

t curves are predicted inverse Compton fluxes for 3 dijfer- 
pulsar wind shock. Hence, measurement JT >x 

., enct assumed B fields. 

of the TeV spectrum over a wider energy 

range may probe spatial variations in the 

nebular magnetic field (see Paper IE). 
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ABSTRACT 

We present results of the analysis of 35.4 hours of observation of the Crab Nebula taken with the 
Whipple Observatory high resolution gamma-ray imaging telescope between October 1995 and 
March 1996. Although we see a steady flux of gamma rays from the nebula, we find no statistically 
significant evidence in the data for emission pulsed at the period of the pulsar. Analysis of 10 hours 
of data obtained for Geminga during the 1996-97 season is in progress and will also be presented at 
the conference. 

INTRODUCTION 

The pulsar in the Crab Nebula and the surrounding nebula is one of the most closely studied 
astrophysical sources. It has been studied extensively in regions of the electromagnetic spectrum as 
diverse as radio waves and TeV gamma rays. In most regions of the spectrum, the 33 ms pulsar 
period is clearly visible when emission from the Crab Nebula is studied. The study of the pulsed 
emission, with its double peaked profile, in different energy ranges is of central importance in the 
development of pulsar models (Eikenberry and Fazio, 1997). The EGRET experiment has shown 
that there is pulsed gamma-ray emission at GeV energies (Ramanamurthy et al. 1995). The Crab 
Nebula is firmly established as a persistent dc source of TeV gamma radiation (Weekes et al. 1989; 
Vacanti et al. 1991; Akerlof et al. 1990; Baillon et al. 1991; Goret et al. 1993). However, these 
observations show no modulation of the signal at the period of the pulsar. In contrast to these reports 
of steady emission from the source, other groups have reported TeV emission modulated at the 33 ms 
period of the Crab pulsar. Some of these reports describe episodic activity (Gibson et al. 1982; Bhat et 
al. 1986; Acharya et al. 1992) while weak pulsed emission persistent over a year has also been 
reported (Do wth waite et al. 1984). 

OBSERVATIONS AND SELECTION OF GAMMA-LIKE EVENTS 

The database analyzed comprised 77 on-source scans with a total duration of 35.4 hours. 55 of these 
scans were taken in an ON/OFF mode under good sky conditions and were used in both the dc analysis 
and the periodicity analysis. The remaining 22 scans, taken in a tracking mode where there were no 
off-source control observations made, were used in the periodicity analysis only. 




All observations were made using the Whipple Observatory 109 element high resolution 
camera at the focus of the 10m optical reflector (Cawley et al. 1990). The camera recorded a 
digitised image of the Cherenkov light pool of each detected shower along with the arrival time of the 
shower which was recorded with a resolution of 0.25 ms. After the standard Whipple methodology 
for flat fielding, gains, normalisation and noise rejection was applied (Fegan et al. 1995), gamma- 
ray-like showers were selected on the- basis of their shape and orientation using the gamma-ray 
selection procedures discussed below. The selection criteria used in the Supercuts-95 and small 
events analyses used are listed in Table 1 while those used in the Extended Supercuts analysis are 
described in detail in Mohanty et al. (1997). 


Table 1. Summary of Event Selection Criteria 


Selection type 

Supercuts-95 

Small events cuts 

Pre-selection 

size > 400 d.c. 

size < 400 d.c. 


maxi > 100 d.c. 

maxi > 40 d.c. 


max2 > 80 d.c. 

max2 > 40 d.c. 


frac3 < 0.975 

frac3 < 0.975 

Distance cut 

0.51°- 1.10° 

0.51°- 1.10° 

Shape cut 

0. 1 6° < length < 0.30° 

length/size < 7.41 xlO' 4 °/d.c. 

0.073° < width <0.15° 

0. 159° < length < 0.269° 
0.080° < width < 0. 142° 

Orientation cut 

alpha < 1 5° 

alpha < 14.7° 

Energy 

> -300 GeV 

-200 - -300 GeV 

Reference 

Catanese et al. (1996) 

Moriarty et al. (1997) 


The standard gamma-ray selection procedure for this data, known as Supercuts-95, was 
optimised on data taken on the Crab Nebula in the 1994-95 season to give the optimum ratio of 
gamma-ray showers selected to statistical fluctuation in the numbers of hadronic images selected 
(Catanese et al. 1996). This procedure incorporates a pre-selection filtering which raises the 
effective threshold energy of the telescope to -300 GeV. 

While the introduction of such a pre-selection filter is desirable from the point of view of 
optimising overall signal to noise ratio, it automatically rejects all small showers. In the context of a 
search for periodic emission from the Crab pulsar which is a source of gamma rays at GeV energies 
(Ramanamurthy et al. 1995) this is clearly undesirable. Accordingly, a separate analysis technique, 
developed to study events in the -200 GeV to -300 GeV region (Moriarty et al. 1997), was used to 
study events with size less than 400 d.c. (digital counts). The most notable difference between this 
small event selection process and Supercuts-95 is the introduction of a length over size cut to reject 
arcs of Cherenkov light rings arising from single muons which are prominent in smaller events. 

Supercuts-95 is optimised to maximise the overall signal to noise ratio but in doing so it 
rejects many of the larger gamma-ray events. A third selection process known as Extended 
Supercuts (Mohanty et al. 1997) was used to help improve the sensitivity of the system to these 
higher energy events. This procedure is quite similar to Supercuts-95 but scales the various cuts 
used with the shower size and passes more of the gamma-ray-like events at higher energies. Since 
the purpose of using Extended Supercuts here was specifically to improve the sensitivity at higher 
energies, a non-standard pre-selection filter of size greater than 400 digital counts was applied. Use 
of Extended Supercuts is appropriate in the context of a Crab pulsar periodicity search since the 
published detections of pulsed emission by Cherenkov telescopes are at energies greater than 1 TeV. 













ANALYSIS AND RESULTS 
DC analysis 

The 55 ON/OFF pairs were analyzed to establish the presence of a steady unpulsed TeV gamma-ray 
signal from the direction of the source over the period of observation. The excess obtained in the 
ON data over the OFF data after the selection criteria were applied was spread through all ON/OFF 
pairs. The results of this analysis are summarised in Table 2. The on-source observing time was 
1521 sidereal minutes and the detector collection area was 3.5xl0 8 cm 2 . 


Table 2. Results of dc analysis 



Raw data 

Supercuts-95 

Small events 

Extended 

Supercuts 

Non 

1321593 

3923 

2374 

7414 

Noff 

1311351 

1607 

1555 

4825 

Non ■ Noff 

10242 

2316 

819 

2589 

Effect 

6.3a 

31.1a 

13.1a 

23.4a 

Integral flux (xl0' n cm' 2 s' x ) 


7.3 


8.1 

Energy 


> -300 GeV 

-200-300 GeV 

>-300 GeV 


Barvcentering and periodicity analysis 

All event arrival times were corrected to the Solar System barycentre using the JPL DE200 
ephemeris (Standish, 1982). The absolute Crab pulsar phase of each event relative to the maximum 
of the radio pulse was then computed using the 
Jodrell Bank Crab pulsar radio ephemeris 
(Lyne & Pritchard, 1997). The precision of the 
barycentering process and the fidelity of GPS 
times appended to events were checked using 
optical observations of the Crab pulsar made on 
December 2 1996 using an infrared photometer 
on the Multiple Mirror Telescope. The output 
pulses from the photometer were fed to the 
Cherenkov telescope data acquisition system 
where a scaler counted the number of 
photometer pulses in 1 ms intervals. The Crab 
pulsar phase for each of the times attached to 
the 1 ms counts was calculated as outlined 
above. Figure 1 shows the phaseogram Pulsar phase 

obtained when the distribution of counts are Fig. 1. Optical light curve obtained after event 
sorted in a 50 bin phaseogram. The clear times are corrected to the Solar System barycentre. 
detection of the pulsar signal, with the main peak in phase with the radio pulse, demonstrates the 
validity of the GPS clock system, data acquisition system and barycentring procedures. 

The different gamma-ray selection criteria discussed above were applied separately to each of 
the 77 on-source scans. Separate 50 bin phaseograms were constructed for each type of cut for each 
data file and for various combinations of data file for each cut (nightly, dark run, full season). Each 
phaseogram was tested for non-uniformity using a % 2 test. There was no statistically significant 
evidence of pulsar activity found in any of the light curves. The lowest chance probability obtained for 
any of the 50 bin phaseograms was 0.006. This was for an Extended Supercuts selection from a single 
run. Given the number of light curves (77 scans, 24 nights, 6 dark runs, full season) and the number of 
different cuts, this value is not statistically significant. Analyses using 25 bin and 10 bin phaseograms 
also yield no significant excesses. 





To calculate upper limits for pulsed emission some assumption has to be made as to the 
pulsar duty cycle. If VHE emission occurs only at phases within the main peak and the EGRET 
main pulse width is 0.1 in phase, the duty cycle is 10%. A 10% duty cycle centred on phase 0.0 was 
assumed and the method of Helene (1983) was used to calculate 99.9% confidence level upper 
limits for persistent pulsed emission over the season of observation. The 99.9% upper limits for 
Supercuts-95 and Extended Supercuts were 76.5 and 95.6 events respectively. The on-source 
observing time was 2125 sidereal minutes. This gave upper limits of 1.7x10' cm' s' for 
Supercuts-95 and 2.1xl0' 12 cm'V for Extended Supercuts with E > -300 GeV in both cases. Thus 
the 99.9% confidence limits based on a 10% duty cycle give a pulsed fraction of less than -2.5% of 
the dc flux. 

DISCUSSION 

The upper limits for pulsed emission presented in this paper improve on the results of an analysis of 
data taken during the 1993-94 observing season and reported at the 24 th ICRC by Gillanders et al. 
(1995). Incorrect 3a upper limits for pulsed emission inadvertently appeared in the text of that paper. 
The correct values presented at the conference for two separate datasets were 3.8xl0' 12 cm'V (E > 
-400 GeV) and 4. lxlO' 12 cm'V (E > -250 GeV). The 3a upper limit on the pulsed fraction for the 
1993-94 observations was -4% of the detected dc flux. The low upper limits obtained for TeV 
pulsed fraction in both these analyses are consistent with results from EGRET (Ramanamurthy et al. 
1995) where the pulsed fraction falls off with increasing energy at energies greater than 1 GeV. 
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ABSTRACT 

PSR B 195 1+32 is ay-ray pulsar detected by the Energetic Gamma Ray Experiment Telescope (EGRET) 
and identified with the 39.5 ms radio pulsar in the supernova remnant CTB 80. The EGRET data 
shows no evidence for a spectral turnover. Here we report on the first observations of PSR B 195 1 +32 
beyond 30 GeV. The observations were carried out with the 10m y-ray telescope at the Whipple Ob- 
servatory on Mt. Hopkins, Arizona. In 8.1 hours of observation we find no evidence for steady or 
periodic emission from PSR B 1951+32 above ~ 260 GeV. FLux upper limits are derived and com- 
pared with model extrapolations from lower energies and the predictions of emission models. 

INTRODUCTION 

The pursuit of Very High Energy (VHE) astrophysics has resulted in the discovery of five sources, 
of which three are associated with young spin-powered pulsars. VHE emission has been detected 
from the direction of the Crab Nebula (Vacant! et.al., 1991), the Vela pulsar (Takanori 1996) and PSR 
B 1706-44 (Kifune et al., 1995) but no evidence has been found for periodic emission at these energies 
in these experiments. 

PSR B 195 1+32 has been detected as a pulsating X-ray source (Safi-Harb et al., 1995) and as 
a high energy y-ray pulsar at E > 100 MeV at the radio period (Ramanamurthy et al., 1995). It can 
be inferred from the five pulsars seen in the MeV to GeV y-ray region that longer period or older 
( ~ 10 5 years ) pulsars have a greater fraction of spin down energy emitted as high energy y-rays. 
The best fit outer gap model of Zhang and Cheng (1997) suggests that PSR B 195 1+32 should emit 
detectable levels of TeV y-rays (Figure 2). The multi wavelength spectrum of PSR B 195 1+32 (Figure 
lb) indicates a maximum power per decade at energies consistent with a few GeV and still rising at 
10 GeV. These factors make PSR B 195 1+32 a good candidate for observations with the ACT above 
100 GeV. 

OBSERVATIONS 

The observations of PSR B 195 1+32 reported here were acquired with the 10 m reflector located at the 
Whipple Observatory on Mt. Hopkins in Arizona. A total of 14 Tracking runs and 4 On /Off pairs 
taken between 13th May, 1996 and 17th July, 1996 constitute the database for all subsequent discus- 
sion. The total On source observing time is 8.1 hrs. The radio position (J2000) of PSR B 195 1+32 (a 
_ 52 m 53 25 j ( 5 = 32° 52' 40.9") was assumed for the subsequent timing analysis. 



Table 1 : Pulsar Parameters 


PSR 

P 

msec 

P 

10 -15 ss -1 

Distance 

kpc 

logioB 

Gauss 

Logio£ 

ergs/s 

B 195 1+32 

39.53 

5.8494 

2.5 

11.69 

36.57 


ANALYSIS AND RESULTS 
Standard Analysis 

The event selection criteria are collectively called Supercuts95 and a detailed description can be found 
elsewhere (Catanese et al., 1995). Supercuts95 raises the effective energy threshold of the detector 
with its software trigger and size cuts. PSR B 195 1+32 appears to have a steep spectrum at EGRET 
energies and since the pulsar spectrum is expected to cut off, it behooves us to reduce the threshold 
of our analysis to search for a lower energy signal. The dominant background at lower energies is 
due to muons whose images appear in the camera as arcs and can be discriminated by a cut on their 
large length/ size values. Hence the selection criteria used Supercuts95 on images with sizes larger 
than 400 p.e. and Smallcuts (Table 2) for images with sizes less than 400 p.e. No steady emission is 

Table 2: Parameter ranges for selecting y-ray images 


Parameter 

Supercuts95 

Smallcuts 

length 

0°16 - 0°30 

unchanged 

width 

0°073 - 0°15 

unchanged 

distance 

0°51 - 1°1 

unchanged 

alpha 

< 15° 

unchanged 

maxi 

> 100 p.e. 

45 p.e. - 100 p.e. 

max2 

> 80 p.e. 

45 p.e. - 80 p.e. 

size 

> 400 p.e. 

0 - 400 p.e. 

length/size 

not used 

< 7.5 x 10" 4 °/p.e. 


Table 3: Selected Events for Steady Emission analysis 


Selection 

Source Events 
a < 15° 

Background Events 
a < 15° 

Excess 

Significance 

Supercuts95 

292 

254 

38 

1.16a 

Smallcuts 

618 

672 

-54 

-1.10a 

Supercuts95 + Smallcuts 

910 

926 

-16 

-0.24a 


detected from PSR B 1 95 1 +32 and 3a flux upper limits are displayed in Table 4. The effective area for 
Supercuts95 , that was used to calculate the upper limit, was taken as A e g ~ 3.5 x 10 8 cm 2 : the same 
area was used for the dataset that resulted from a combination of Supercuts95 and Smallcuts although 
here there is more systematic uncertainty. The energy threshold was obtained from simulations and 
extrapolating the Crab Nebula y-ray rate for each set of cuts used assuming a spectrum ~ E -2 4 . 

Periodic Analysis 

The arrival times of the Cerenkov events were registered by a GPS clock with an absolute resolution 
of 250 /isec. An oscillator calibrated by GPS second marks was used to interpolate to a resolution of 






0. 1 /isec. After an oscillator drift correction was applied, all arrival times were transformed to the solar 
system barycenter and folded to produce the phases, 4>y, of the events modulo the pulse period. The 
ephemeris frequency parameters used were v= 25.2963719901267 s -1 and v=-3.73940x 10 12 s 2 , 
at the epoch to=JD 2450177.5. This frequency was extrapolated 72 days to obtain a timing solution 
relevant to the epoch of observation. The datasets, however, were taken within the validity interval of 
the above ephemeris. 

To check the Whipple Observatory timing systems an optical observation of the Crab pulsar 
was undertaken on December 2nd (UT) 1996 using the 10m reflector. The phase analysis of the event 
arrival times yielded a clear detection of optical Crab pulsar signal which is in phase with the radio 
pulse and demonstrates the validity of the timing, data acquisition and software in the presence of a 
pulsed signal. No evidence of pulsed emission from PSR B 195 1+32 at the radio period exists. To 
calculate a pulsed flux upper limit we assumed the same pulse profile as seen at EGRET energies, i.e. 
with the phase range for the main pulse and secondary pulse as 0. 12 - 0.22 and 0.48 - 0.74 respectively 
(Ramanamurthy et al., 1995). 

Table 4: Integral Flux Upper limits 



Steady Emission 

Periodic Emission 

Threshold 


cm -2 s -1 ) 

(cm -2 s _1 ) 

(GeV) 

Supercuts95 

0.97 x 10" 11 

3.7 x 10~ 12 

> 370 

Supercuts95 + Smallcuts 

1.95x 10' 11 

6.7 x 10' 12 

> 260 


DISCUSSION 

PSR B 195 1+32 is surrounded by a compact nebula 
which gives a plerionic nature to the supernova rem- 
nant, CTB80. X-ray plerions are good candidates for 
VHE emission since the electrons responsible for neb- 
ular synchrotron X-rays should also create VHE y-rays 
via the inverse Compton (IC) process. It is expected 
that for plerions, such as that associated with PSR 
B 195 1+32 where the density of nebular synchrotron 
photons is too low for SSC to take place, detectable 
VHE emission should be produced by the IC scatter- 
ing of the 2.7K cosmic microwave background by the 
same electrons radiating synchrotron X-ray photons. 
Interpreting the unpulsed X-ray emission form CTB80 
as the synchrotron emission from a plerion, the esti- 
mated IC flux > 1 TeV is 6.6 x 1CT 13 TeV/cm 2 /srteV 
(De Jageret. al., 1995). This represents the lower limit 
on the IC flux since there can be other sources of soft 
photons in addition to the microwave background. 

To model the pulsed high energy spectrum, a 
function of the form 

dNy/dE = KE~ r e 



Fig. 1: The pulsed energy spectrum of PSR 
B 195 1+32. The Whipple limit is indicated 
as a filled square at 370 GeV. (See text for 
details). 

i-E/Eo) (J) 


was used where E is the photon energy, T is the photon spectral index and Eo is the cut off energy. 
The pulsed upper limit reported here is two orders of magnitude lower than the extrapolated EGRET 
power law. Equation (2) was used to extrapolate the EGRET spectrum to VHE energies constrained 





by the TeV upper limit reported here and indicates a cut off energy of Eq < 75 GeV for pulsed 
emission (Figure la). 

The strength of the cut off provides a good dis- 
criminant between the various pulsar emission mod- ^ 10 
els. The status of current observations and the derived * 5 
cutoff discussed above indicates that the cutoff is be- T g 
yond 10 GeV. In polar cap models this would indicate j 
a sharp cutoff since the pair production optical depth 5 " 5 
increases exponentially with photon energy (Harding -to 
1997). However, it is not possible to constrain the _, s 
shape of the cut off with the non detection of pulsed 
TeV flux reported here. The most relavant compari- 
son of the Whipple upper limit with emission models 
is the outer gap model of Zhang and Cheng (see Fig- 
ure 2). This model includes the effect of geometry in _ 35 
the treatment of pulsed emission via a parameter a = 

.i . . -40 

r/rj,, the radial distance to the synchrotron emitting re- 
gion near the outer gap, r, as a function of the light 
cylinder radius r L . Our pulsed upper limits are consis- Fig 2 : Predicted pulsed y-ray flux of PSR 
tent with the outer gap model if a > 0.6 implying an B1951+32from the Zhang and Cheng outer- 
emission region far out in the magnetosphere. gap model. The solid, dot-dash and dashed 

The result reported here is the first observation curves correS p 0 nd to a =0.5, 0.6, 0.7 respec- 
of PSR. B1951+32 beyond 30 GeV. PSR B1951+32 tively. (See text for details). 
exhibits very similar spectral behavior and morpho- 
logical features, such as an associated synchrotron nebula, to PSR B1706-44 (Finley et al., 1997). If 
these factors are any indication of similar emission mechanisms in pulsars then the lack of unpulsed 
emission from PSR B 195 1+32 is puzzling considering that PSR B 1706-44 was detected as a VHE 
source of unpulsed emission > 1 TeV (Kifune et al., 1995). Lack of pulsed emission indicates that the 
processes producing pulsed high energy photons over two decades of energy in the EGRET energy 
range somehow become ineffective over a decade of energy to result in a lack of VHE y-rays. The 
low magnetic field of PSR B 195 1+32 relative to the average pulsar field implies that attenuation of 
y-rays by magnetic absorption is not a likely explanation for the non-detection. 
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ABSTRACT 

If supernova remnants (SNRs) are the sites of cosmic-ray acceleration, the associated nuclear interac- 
tions should result in observable fluxes of TeV gamma-rays from the nearest SNRs. Measurements of 
the gamma-ray flux from six nearby, radio-bright, SNRs have been made with the Whipple Observa- 
tory gamma-ray telescope. No significant emission has been detected and upper limits on the >300 
GeV flux are reported. Three of these SNRs (IC443, gamma-Cygni and W44) are spatially coincident 
with low latitude unidentified sources detected with EGRET. These upper limits weaken the case for 
the simplest models of shock acceleration and energy dependent propagation. 


INTRODUCTION 

It is generally believed that cosmic rays with energies less than ~ 100 TeV originate in the galaxy and 
are accelerated in shock waves in shell-type SNRs. This hypothesis is supported by several strong 
arguments. First, supernova blast shocks are one of the few galactic sites capable of sustaining the 
galactic cosmic ray population against loss by escape, nuclear interactions and ionization energy loss 
assuming a SN rate of about 1 per 30 years and a 10% efficiency for converting the mechanical 
energy into relativistic particles. Second, models of diffuse shock acceleration provide a plausible 
mechanism for efficiently converting this explosion energy into accelerated particles with energies 
~ 10 14 - 10 15 eV and naturally give a power-law spectrum similar to that inferred from the cosmic 
ray data after correcting for energy dependent propagation effects. Finally, observations of non- 
thermal X-ray emission in SN1006 (Koyama, et al., 1995) and IC443 (Keohane, et al., 1997) suggest 
the presence of electrons accelerated to ~ 100 TeV and ~ 10 TeV respectively. 

If SNRs are sites for cosmic ray production, there will be interactions between the acceler- 
ated particles and the local swept-up interstellar matter. Drury, Aharonian and Volk (1994) (DAV) 
and Naito and Takahara (1994) have calculated the expected gamma-ray flux from secondary pion 
production using the model of diffusive shock acceleration. The expected intensity (DAV) is 


/ E \ / 9 Esn 

F(>E) = 9x10 1TeV J y 10 51 erg 


lkpcy Vlcm 


-3) 


cm 


( 1 ) 


where E is the photon energy, 0 is the efficiency for converting the supernova explosion energy, E S n , 
into accelerated particles, d is the distance to the SNR and n is the density of the local ISM. 



OBSERVATIONS 

We report on the results of observations of six nearby SNR (W44, W51, gamma-Cygni, W63, Ty- 
cho and IC443) by the Whipple Observatory’s high energy gamma-ray telescope situated on Mount 
Hopkins in southern Arizona. The telescope (Cawley et al., 1990) employs a 10 m diameter optical 
reflector to focus Cerenkov light from air showers onto an array of 109 photmultipliers covering a 3 
degree field of view. By making use of distinctive differences in the lateral distribution of gamma-ray 
induced showers and hadronic induced showers, images can be selected as gamma-ray like based 
on their angular spread. The determination of the incident direction of the selected gamma-ray like 
events is accomplished by making use of the orientation, elongation and asymmetry of the image. 
Monte Carlo studies have shown that gamma-ray images are a) aligned towards their source position 
on the sky b) elongated in proportion to their impact parameter on the ground and c) have a cometary 
shape with their light distribution skewed towards their point of origin in the image plane. Results on 
the Crab Nebula indicate that the angular resolution function for the telescope using this technique is 
a Gaussian with a standard deviation of 0. 13 degrees. (Lessard et al., 1997). A combination of Monte 
Carlo simulations and results on the Crab Nebula indicate that the energy threshold of the technique 
is 300 GeV and the effective collection area for a point source located at the center of the field of view 
is 2. 1 x 10 8 cm 2 and is reduced for offset sources (Lessard et al., 1997). 


The analysis of data from extended sources involves 
binning the event arrival directions. We define the source 
region for the SNR by a circular aperture which matches the 
maximum extent of the radio shell (Green, 1995) plus twice 
the width of the angular resolution function to account for 
the smearing of the edge of the remnant. The number of 
gamma-ray candidate events is obtained by subtracting the 
number of events in the OFF-source observations from the 
number of events in the ON-source observations. 

RESULTS 
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The observations were made over three observing seasons, 
from 1993 to 1996. Two dimensional images of the excess 
events for the Crab Nebula (which was deliberately offset 
from the center of the camera, to demonstrate the verac- 
ity of the technique) and the SNR W51 are shown in Fig- 
ure 1. In each frame, the statistical significance is displayed 
in grayscale. The black contours are from the 4850 MHz ra- 
dio survey by Condon et al., 1994, showing the extent of the 
radio shell. The circle shows the circular aperture used to de- 
rive the excess counts from the entire remnant (see Table 1). 

Of the six remnants observed, W 51 showed the great- 
est excess at an offset location which is spatially consistent 
with hard x-ray emission and peak radio intensity. More data 
are required to determine if this excess is signicant. No sig- 
nificant excess has been recorded for the other remnants and 
99.9% confidence upper limits on the flux have been calcu- 
lated (see Table 1). The upper limit assumes uniform emis- 
sion from the remnant in the absence of a priori knowledge 
of a more defined emission region. 

DISCUSSION 

In Figure 2 the Whipple upper limits and EGRET data (Es- 
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Fig. 1: Observations of the Crab Neb- 
ula and W51. We have applied a 
boxcar smoothing method which max- 
imizes the point source sensitivity. The 
statistical significance of the excess se- 
lected events is shown in grayscale. 
The overlayed contours are propor- 
tional to the 4850 MHz radio intensity. 


posito et al., 1996, for gamma-Cygni and IC443, Fierro, 1995, for W44, and Thompson et al., 1995 
for the remaining) are compared with an E~ 2 * extrapolation of the EGRET data using the contri- 
bution to the gamma-ray spectrum from secondary pion decay as derived by Buckley et al. (1997) 
using the model of DAV. The upper dotted curve assumes a source spectrum of E~ 2A and a reasonable 
maximum value of the product Es/yQ j df used in the model. The lower dotted curve assumes a source 
spectrum of E~ 2A and a reasonable minimum value of the product E$nQ / d . 

We interpret our results in the 
context of two hypotheses, (1) that the r uf* 

EGRET data gives evidence for accel- •’ 10 _, 
eration of cosmic ray nuclei in SNR s )0 -» 
and that the observed gamma-ray emis- 1 10 - s 
sion comes not from primary electrons 
but from nuclear interactions of cosmic l, 10 . , 2 

w 1 0 

rays with ambient material or (2) that - |Q -u 
the EGRET flux is produced by some 10 - u 
other mechanism. 10 

Under the assumption that the 



contribution from electron bremsstra- 
hlung and inverse Compton (IC) scat- 
tering are neglible, it is reasonable to 
compare the high energy gamma-ray 
upper limits to an extrapolation of the 
integral EGRET fluxes using the model 
by DAV. In the case of gamma-Cygni, 
IC443 and W44 the Whipple upper lim- 
its lie a factor of ~ 25, 10 and 10 re- 
spectively below the extrapolation and 
require either a spectral break or a 
source spectrum steeper than E" 2 5 for 
gamma-Cygni and £ -2 4 for IC443. 

Another plausible explanation for 
the results is that the EGRET flux is 
produced by high energy electrons ac- 
celerated in the vicinity of pulsars. If 
this is the case, then the Whipple up- 
per limits must be compared with the 
a priori model predictions. There is 



Energy (CeV) 



Energy (CeV) 


enough uncertainty in the parameters of 2.- Whipple upper limits shown along with EGRET in- 
the SNR that the upper limits are not in t e g r al fluxes, and integral spectra. These are compared to 
strong conflict with these predictions, extrapolations from the EGRET integral data points (solid 
but it is still strange that in these ob- curve s) i as well as a conservative estimate of the allow- 
jects which show strong evidence for a ^ e ran g e of fluxes from the model of DAV ( dotted curves). 
interactions with molecular clouds (cor- ^ SQ s ^ own are CASA-MIA upper limits from Borione et al., 
responding to the upper dotted curve) 7995 CYGNUS upper limits from Allen et al., 1995, and the 
in no case is there an observable TeV A jpQBlC upper limit from Prosch et al, 1996. 
gamma-ray flux. Evidence of an X-ray 

point source embedded in gamma-Cygni (Brazier et al., 1996) and IC443 (Keohane et al., 1997) and 
the observation of a pulsar, B 1853+01, in W44 (Wolszczan, et al. 1991), all provide support to a 


pulsar origin for the EGRET flux. 









Table 1 : Results of Observations. 


Source 

Pointing 

Aperture 

Name 

a, 5 

Radius 


(1950) 

(deg) 

W44 

18:53:29,01:14:57 

0.55 - 

W51 

19:21:30, 14:00:00 

0.68 

y-Cygni 

20:18:59,40:15:17 

0.76 

W63 

20:17:15,45:24:36 

1.05 

Tycho 

00:22:28, 63:52:11 

0.29 

IC443 

06:14:00, 22:30:00 

0.64 


ON 

OFF 

Total 

Upper 

Source 

Source 

Time 

Limit x 10' 

Counts 

Counts 

(min) 

(cm -2 s -1 ) 

450 

426 

360.1 

3.0 

619 

559 

468.0 

3.6 

1040 

1104 

560.0 

2.2 

452 

501 

140.0 

6.4 

315 

302 

867.2 

0.8 

1565 

1522 

1076.7 

2.1 
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ABSTRACT 

The ~70 unidentified sources of the EGRET sky survey may be one of its most important legacies. 
The identification of these sources at other wavelengths is critical to understanding their nature. Many 
have flat spectra out to 10 GeV which, if extrapolated to TeV energies, would be easily detectable 
relative to the steeply falling diffuse background. The Whipple Observatory y-ray telescope has been 
used to observe a number of these which were selected based on their position, intensity and spectrum 
and in some cases based on a possible association with a supernova remnant or pulsar. No significant 
emission has been detected from these sources, and upper limits are given. 

INTRODUCTION 

Despite extensive searches for counterparts to the ~30 low-lattitude unidentified EGRET sources, 
the nature of these objects is still largely unknown. Kaaret & Cottam (1996) have suggested that the 
low latitude unidentified sources show a correlation with OB associations, the sites of massive star 
formation. Since nearly half of all supemovae occur as the core collapse of young massive stars which 
explode into star formation regions (e.g., Huang & Thaddeus 1986) a correlation with the positions 
of pulsars and with supernova remnants also follows. Kaaret & Cottam argue that pulsars emit y- 
rays over a significantly longer lifetime than SNR, so that the number of y-ray pulsars is expected 
to be significantly larger than the number of y-ray SNR. However, Esposito et al. (1996), Stumer & 
Dermer (1994) and Stumer, Dermer and Mattox (1996) have presented evidence for associations of a 
number of these objects with SNRs (y-Cygni, IC443, W44, Monoceros) for which there is no pulsar 
within the 95% confidence error contour (Stumer, Dermer & Mattox 1996). Attempts to detect radio 
pulsars in the error boxes of EGRET unidentified sources have been unsuccessful (e.g., Nice & Sayer 
1997) and provide some constraints on models in which all of the Galactic unidentified sources are 
pulsars. Since EGRET generally lacks the spatial resolution to distinguish the point-like emission 
from pulsars and AGNs from the extended emission expected to arise in the vicinity of supernova 
shells, variability has been used to distinguish compact sources. Dramatic transient sources such as 
the enigmatic 2CG 135+1 and newly discovered GRO J 1838+04 are difficult to interpret as either 
arising from AGNs or from pulsars, and are possibly representatives of a new class of Galactic y-ray 
source distinct from isolated pulsars (e.g., Tavani et al. 1997). 



Detection of these objects at high energies using ground based y-ray detectors would aid in the 
identification of these objects in two important ways. First, the contribution from the diffuse y-ray 
background falls as a steeper power of energy (~ E~ 21 ) than the source spectrum (~ E~ 2 ) for many 
of these objects, implying a smaller effect from uncertainties in the diffuse background model in 
determining the position, flux and spectra of these sources. Second, the 0.13° angular resolution of 
the Whipple 10m y-ray telescope (Lessard and Buckley 1997) provides the ability to resolve extended 
sources such as SNRs and offers the potential to narrow the error box for bright sources. 

OBSERVATIONS AND ANALYSIS 

The high-energy y-ray telescope (Cawley et al. 1990) at the Whipple Observatory employs a 10 m 
diameter optical reflector to image Cerenkov light from air showers onto an array of 109 fast pho- 
tomultipliers covering a 3° field of view (FOV). By making use of the distinctive differences in the 
angular distribution of light and orientation of the shower images a y-ray signal can be extracted from 
the large background of hadronic showers. 

Data are generally taken in a differential mode where each 28 min ON-source run is followed 
by a 28 min OFF-source control run which is offset in right ascension to ensure that the same range 
of azimuth and zenith angles are sampled. While this cancels the zenith angle dependence of the 
cosmic-ray rate as well as other systematic effects in the camera, differences in sky brightness between 
the ON-source and OFF-source regions can lead to biases. For some galactic plane sources, such 
differences are substantial due to either diffuse emission from the galactic plane or bright stars. 

Systematic effects arising from such brightness differences can be largely canceled by the pro- 
cedure of software padding (Cawley et al. 1983). This procedure consists of adding noise to all pixels 
of each event so that matching PMTs in ON-source and OFF-source runs have identical noise pulse 
height spectra. Only PMT signals which exceed some multiple of this noise level are included in the 
subsequent analysis of the shower images (Punch et al. 1993). 

The technique used to generate the two-dimensional plots and upper limits is a simple extension 
of the standard Whipple data analysis (Reynolds et al. 1993) and is described in more detail in Lessard 
and Buckley (1997). After initial processing of the shower images including pedestal subtraction, 
gain correction and image cleaning (e.g.. Punch et al. 1993) individual Cerenkov shower images are 
subjected to a moment analysis to determine a set of parameters that characterize the roughly elliptical 
images. Each point on a two dimensional grid covering the 3° FOV is considered as the potential 
source position. For each event, the RMS width and length, centroid position, orientation, ellipticity 
and the skew of the shower image are calculated about this point of origin and tested for consistency 
with the parameter values expected for a y-ray event coming from the corresponding direction in the 
sky. For each grid point the number of candidate events ON-source and OFF-source are calculated, 
and the significance of the excess, S\, is derived using the likelihood ratio method of Li & Ma (1983). 
In the two-dimensional plots presented in Figures 1, the gray-scale indicates the number of excess 
counts (candidate y-rays) consistent with each grid point and the contours shown correspond to the 
statistic S\ in steps of 1. Note that due to the large number of trials associated with the 30x30 grid, 
approximately one S \= 3 excess is expected for each two-dimensional plot. 

While it is desirable to have one control (OFF-source) run for each run ON-source, for some 
of the data presented here the number of exposures taken OFF-source is less than the number of 
ON-source runs. In this case, the background level is determined by normalizing the OFF-source 
data to the ON-source data in a 0.25° band around the perimeter of the field of view. The resulting 
normalization factor a enters into the calculation of the significance and the upper limit following the 
procedure of Li and Ma (1983). For the sources J0542+26, J0635+0521, and J 1825- 1307, little or no 
OFF-source data was taken and a background template was formed using contemporaneous control 
data taken for other sources. This results in an additional systematic error for these sources. 

In calculating upper limits, we are testing the hypothesis that the emission is coming from a 
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point source within the EGRET error box. Upper limits calculated for the extended regions corre- 
sponding to the IC443, W44 and y-Cygni SNRs (for 2EG J06 18+2234, 2EG J 1857+01 18, and 2EG 
J2020+4026 respectively) are reported elsewhere (Buckley et al. 1997). 99.9% confidence-level upper 
limits are calculated for each grid point lying within the EGRET 95% confidence interval using the 
method of Helene (1983) and accounting for the declining y-ray detection efficiency away from the 
camera center (Lessard & Buckley 1997). The maximum upper limits for each error box are shown 


in Table 1. 


RESULTS AND CONCLUSIONS 

Table 1 shows preliminary Whipple upper limits for a number of unidentified sources together with 
the extrapolated EGRET flux derived from the EGRET spectrum. Fluxes and spectral indices are from 
Fierro et al. (1996). In addition to sources from the first (Fichtel et al. 1994) or second EGRET cata- 
logs (Thompson et al. 1995), we also include the source J0749+17 from the initial list of unidentified 
sources by Hartman et al. (1992). This object was not included in the first EGRET catalog because of 
its low significance « 4a), but is of interest since it prompted the radio pulsar search by Lundgren, 
Zepka and Cordes (1995) that led to the discovery of the binary millisecond pulsar PSR 0751+1807. 
Also on our list is the source J0542+26 which was on the list of high confidence unidentified sources 
in the first but not the second EGRET catalog. This source has a 158 arcmin error radius which could 
not be shown in Figure 1 . This source is of interest since it is coincident with the position of the old, 
nearby (0.8-1.4kpc, Kundu et al. 1980) SNR S 147 as pointed out by Stumer and Dermer (1994). 

These data were taken over the period December 1994 to May 1997. Two-dimensional plots for 
these sources are shown in Figure 1 excluding results for J06 18+2234, J 1857+01 18 and J2020+4026 
which are shown elsewhere (Buckley et al. 1997). Upper limits are at energies above 400 GeV unless 
otherwise indicated. 2EG 1746-2852 transits at an elevation of <30° resulting in an increase of the 
effective area and energy threshold by a factor of approximately 5.0 compared with observations at the 
zenith (Krennrich et al. 1997). While 2EG J 1746-2852 shows a small (2.5a) excess at the position of 
Sgr A* and within the EGRET error box, this excess is not considered significant given the additional 
systematic errors present for galactic plane sources. 2EG J0241+61 19 shows a similar excess near the 
position of 2CG 135+01 and within the EGRET error box. The excess in J0542+26 lies outside and to 
the south of the radio shell of S 147 and approximately 0.5° away from the X-ray binary 4U0535+262, 
too far to make an association. The other sources show no significant emission within the EGRET 
error boxes. Further deep observations with the GRANITE-m high resolution camera should provide 
better sensitivity given the extended FOV and finer pixelization, and correlations with data taken at 
other wavelengths should improve chances for detecting variable unidentified sources. 
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Fig. I: Plots of two-dimensional distributions of candidate gamma-ray events and contours of S, 
for (a) J024I +61 19 (cross at the position of 2CG1 35+01 ), ( b ) J0542+26, (c) J0545+ 3943. (d) 
J0635+052 1 , (e) J 1746-2852 (cross at Sgr A ') and (f) J 1825+2234 (cross at the position of PSR 
B1 823-1 3). Dotted contours show the elliptical fit to the EGRET 95% confidence intervals. . 
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